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Application Serial No.: 09/739,940 
Examiner My-Chau T. Tran 
Group art Unit: 1639 
Attorney Docket No. 30626-101 

REMARKS/ ARGUMENTS 

The Final Office Action of March 1 1 , 2004 has been carefully reviewed and this 
response addresses the Examiner's concerns. 

Status of the Claims 

Claims 1-3, 6-7, 10-12, 14-19, 66-67, 69-70 and 125 were pending in the application. 
Claims 1 19-124 are withdrawn from consideration. 

In view of the examiner's earlier and present restriction requirements, applicants retains 
the right to present claims 1 19-124 in divisional applications. 

Claims 1-3, 6-7, 10, 12, 14-19, 66-67 and 69-70 and 125 were rejected under 35 USC 
102(b) as being anticipated by Bal et al. (AnaL Chem,, 1994, 66:3423-3430). 

Claims 1-3, 6-7, 10-12, 14-19, 66-67, 69-70 and 125 were rejected under 35 U.S.C. 
§103(a )as being unpatentable over Bal et al. in view of Farmer et al. (J. Mass Spectrom., 1998. 
3:697-704). 

Claim 1 is amended to further describe the Applicants' invention. 

Support in the Specification to Amendments to Independent Claim 1 

The analyzing the sample by matrix-less light desorption/ionization mass spectroscopy, 
after the sample has been applied to the deposited continuous film 1 is supported by p. 23, lines 
19-21 {''The component of the present invention dealing with light (laser) desorption for mass 
spectrometry is a matrix-less laser desorption approach. It uses the deposited semiconductor 
films of Table L p. 29, lines 12-14 (''The proteins, peptides and ammonium citrate used to 
demonstrate matrix-less mass spectroscopy with the three film morphologies of Table 1 were 
obtained from Sigma and the trypsin was frozen, sequencing grade from Promega. p.30, lines 
16-17 and Fig. 4 ("FIG, 4 typifies results from a continuous (void free) film. p. 36, lines 1-2 
("/w light (laser) desorption using the matrix-less approach, rapid thermal heating of non-matrix 
substrates is what releases the molecules from the surface. % 
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Claims 1-3, 6-7, 10, 12, 14-19, 66-67 and 69-70 and 125 were rejected under 35 USC 102(b) as 
being anticipated by Bal et al (Anal Chem., 1994, 66:3423-3430). 

Applicants respectfully traverse these rejections for the reasons presented below. 

In order to better understand the differences between the present invention and the Bal et 
al. paper it should be recognized that in the Applicants' invention the continuous films have 
optical properties and species adsorption properties essentially the same as the optical properties 
and the species adsorption properties of the continuous film bulk material and the sample is 
analyzed by matrix-less light desorption/ionization mass spectroscopy . 

In comparison, as generally pointed out by the examiner, the Bal et al. paper discloses a 
method of analyzing protein by matrix-assisted laser desorption/ionization mass spectrometry 
(MALDI mass spectrometry) (Abstract; pg. 3423, left col, lines 1-6; pg. 3424, left coL, line 45 
to right col., line 17). The method comprises depositing a Nafion film by evaporation onto a 
stainless-steel probe (metal substrate) (pg. 3425, lines 13-27) and applying the sample directly 
onto the Nafion film after the deposition of the film (pg. 3425, left col., line 39 to right coL, line 
2). The Nafion film is a perfluorinated cation exchanger with sulfonic acid groups (carbon 
hydrogen mixture) and possesses a high surface activity and can selectively bind small cations 
and organic cations covalently (species adsorption/analyte adsorption) (pg. 3424, left col., line 45 
to right col., line 17). Nafion in liquid form at a concentration of 5% in a mixture of 10% water 
and alcohols was used to apply to the stainless steel probe. The solution was air-dried to produce 
the film. (pg. 3425, left column, lines 13-17). 

The following summary of the prior art description of the properties of Nafion is 
presented to support that at least one of the Applicants' claimed limitations, specifically, 
depositing a continuous film having optical properties and species adsorption properties 
essentially the same as the optical properties and the species adsorption properties of the 
continuous film bulk material, is not found in the Bal et al. paper. 

"Nafion is initially produced in a salt form with its sulfonic acid groups 
neutralized." (Nafion: Physical and Chemical Properties available at 

http://www.permapure.com/newweb/Nafion%20properties.htm and provided in Appendix A). 
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As understood in the art, in the salt form, Nafion is not chemically active. Nafion can be 
"chemically activated by conversion of the salt form to the acid form." "Once activated, Nafion 
immediately begins to react with its surrounding environment. Moisture is absorbed and 
exchanged with the surroundings. Reactive organic gas components within the surrounding air 
may undergo chemical changes when exposed to the super-acid catalytic activity of Nafion. 
These organic gases may combine to form larger compounds that are liquid or solid in nature. 
Over time, these organic residues build up a deposit on the Nafion. The original color of the 
Nafion changes gradually from translucent to yellow, then brown, then even black." (Nafion: 
Physical and Chemical Properties available at 
http://www.permapure.com/newweb/Nafion%20properties.htm). 

" [The] ionic properties of Nafion were created by adding sulfonic acid groups, a 
chemical with very strong ionic properties, into the bulk polymer matrix, Nafion is a super-acid 
catalyst. The sulfonic acid groups attached to the Teflon backbone within Nafion function as an 
extremely strong proton donor due to the stabilizing effect of the large polymer matrix attached 
to the sulfonic acid. Nafion is very selectively and highly permeable to water. The sulfonic acid 
groups in Nafion have a very high water-of-hydration, so they very efficiently absorb water. 
Interconnections between the sulfonic acid groups lead to very rapid transfer of water through 
the Nafion." (Nafion: Physical and Chemical Properties available at 
http://www.permapure.com/newweb/Nafion%20properties.htm). 

Nafion films are originally amorphous but develop crystalline regions after a short period 
of time (a week). (M. Ludvigsson et al., J. Electrochem. Soc, 2000, 147: 1303-1305, provided in 
Appendix B, see pg. 1304, rt. column, second par.) the crystalline regions in Nafion become 
amorphous when wetted. (Ludvigsson, pg. 1305, left column, first paragraph). Since Nafion 
absorbs water, a Nafion film would have amorphous and crystalline regions depending on the 
water absorption. As the Nafion film absorbs water from the ambient, the film expands both 
vertically and laterally. (P.J. James et al., J. Material Sci., 2000, 35:51 1 1-51 19, provided in 
Appendix C, see pg. 5117, left column, line 5). Both AFM images and X-ray scattering indicate 
that clusters (clusters of ionic aggregation) form in the Nafion film and the clusters vary with 
water content, i.e. humidity (P.J. James, Figs. 3, 5, 6). Furthermore, it has also been shown that 
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the conductivity of the Nafion film varies with water content (Jeffrey B. Mecham, Ph. D. Thesis, 
Department of Chemistry, Virginia Polytechnic Institute and State University, 2001, provided in 
Appendix D, see pg. 84-87). 

The Bal et al. paper is silent as to the optical properties. However, it is inherent that 
materials with different ratios of crystalline to amorphous volumes have different optical 
properties. (See for example, the use of crystalline to amorphous transitions in rewritable optical 
media such as CD-RW. See definition of phase change disk in Webopedia at 
http://www.webopedia.eom/TERM/P/phase_change_disk.html a copy of which is provided in 
Appendix E). It is also inherent that materials with different conductivity have different optical 
properties since the materials interact differently with electromagnetic waves. Furthermore, it is 
also inherent that materials with different topography will scatter electromagnetic waves 
differently and will have different optical properties. 

The Bal et al. paper discloses the use of a Nafion film in MALDI analysis of an analyte. 
In one example the analyte is a protein sample, The protein sample is obtained by electrophoretic 
analysis and HPLC (pg. 3424, right col., line 50 to pg. 3425, left col, line 12; pg. 3425, lines 
30-39). The sample is mixed with a matrix (pg. 3424, right col., lines 37-49) for MALDI. The 
sample is analyzed by laser mass spectrometry (pg. 3424, right col., line 20-36). 

Even further supporting Applicants' position stated above, presented below are details 
related to optical properties and adsorption properties of the Nation film. 
Optical properties 

The bulk material described in the Bal et al. paper is Nafion in liquid form at a 
concentration of 5% in a mixture of 10% water and alcohols. That a solution has different optical 
properties form those of the solute dissolved in the solvent when the solute is considered 
separately or any of the separate components follows from the classical work of Onsager, 
Kirkwood and others (see references in G, van der Zwan, Theories for the dielectric constant, 
available at www.chem.vu.nl/acas/Homepages/zwan/notes.html a copy of which is provided in 
Appendix F). Thus, the Nafion film does not have the same optical properties as the base 
solution at a concentration of 5% in a mixture of 10% water and alcohols. Furthermore, the 
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Nafion film will have different optical properties at different humidities since the Nafion film 
will have amorphous and crystalline regions depending on the water absorption according to 
Ludvigsson et al. while the Nafion solution at a concentration of 5% in a mixture of 10% water 
and alcohols would not have a variable reflectivity due to amorphous and crystalline regions. 

Assume arguendo that the salt form (with its sulfonic acid groups neutralized) is 
considered the bulk material. In the salt form, Nafion is not active. The Nafion film has to be 
active for the film to serve as part of the active substrate for MALDI (Bal et al., pg. 3424, left 
column, last paragraph). The active film will have different optical properties at different 
humidities and organic residues will build up a deposit on the Nafion changing the color (optical 
property) gradually from translucent to yellow, then brown, then even black. Thus, the Nafion 
film does not have the same optical properties as the salt form. 

Adsorption Properties 

The bulk material described in the Bal et al. paper is Nafion in liquid form at a 
concentration of 5% in a mixture of 10% water and alcohols. In a solution the solvent molecules 
interact with the solute molecules in order to form loosely bonded aggregates. (C.W. Keenan, J. 
H. Wood, General College Chemistry, 3rd edition, 1966, p.236, provided in Appendix G). Thus, 
the sulfonic acid groups are not available to act as extremely strong proton donors and can not 
bind cations. In its active form, the sulfonic acid groups act as extremely strong proton donors 
and can bind small and large organic cations. Thus, the species adsorption properties of the 
Nafion film are not essentially the same as the species adsorption properties of the Nafion bulk 
material. 

Assume arguendo that salt form (with its sulfonic acid groups neutraUzed) is considered 
the bulk material. With the sulfonic acid groups neutralized, the sulfonic acid groups are not 
available to act as extremely strong proton donors and can not bind cations. In its active form, 
the sulfonic acid groups act as extremely strong proton donors and can bind small and large 
organic cations. Thus, even under this assumption, the species adsorption properties of the 
Nafion film are not essentially the same as the species adsorption properties of the Nafion bulk 
material. 
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Thus, Bal et al do not teach depositing a continuous film having optical properties and 
species adsorption properties essentially the same as the optical properties and the species 
adsorption properties of the continuous film bulk material. 

Furthermore, the Bal et al. paper discloses a method of analyzing protein by 
matrix-assisted laser desorption/ionization mass spectrometry (MALDI mass spectrometry) (Bal 
et al. Abstract; pg. 3423, left col., hnes 1-6; pg. 3424, left col., line 45 to right col., line 17). As 
disclosed by Bal et al., the success of the MALDI experiments disclosed in the Bal et al. paper 
depends on the choice of the UV absorbing matrix (Bal et al. , p. 3423, rt. column, lines 5-6). In 
contrast, in the Applicants' claimed invention the sample is analyzed by matrix-less light 
desorption/ionization mass spectroscopy. 

Thus, Bal et al. do not teach analyzing the sample by matrix-less light 
desorption/ionization mass spectroscopy, after the sample has been applied to the deposited 
continuous film. 

To anticipate a claim a reference must teach every element of the claim. (MPEP § 
2131). As pointed above, Bal et al do not teach depositing a continuous film having optical 
properties and species adsorption properties essentially the same as the optical properties and the 
species adsorption properties of the continuous film bulk material, and Bal et al also do not teach 
analyzing the sample by matrix-less light desorption/ionization mass spectroscopy, after the 
sample has been applied to the deposited continuous film. 

Applicants respectfully assert that Claim 1 is not anticipated by the Bal et al paper and 
neither are any of the dependent claims. In addition, in view of the specific needs of Bal et al (as 
discussed above) and since Bal et al is lacking at least one patentable feature present in claim 1 
and the dependent claims of the Applicants' invention, a modification of Bal et al under 35 
U.S.C. §103 would also be inappHcable because such modifications are not taught by nor 
obvious under Bal et al, and if incorporated into Bal et al, would render Bal et al. inoperable for 
its intended functions. 
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Claims 1-3, 6-7, 10-12, 14-19, 66-67, 69-70 and 125 were rejected under 35 U.S.C, §103(a )as 
being unpatentable over Bal et aL in view of Farmer et al (J. Mass Spectrom,, 1998, 3:697-704), 

As stated above, Bal et al. do not disclose depositing a continuous film having optical 
properties and species adsorption properties essentially the same as the optical properties and 
the species adsorption properties of the continuous film bulk material and Bal et al. also do not 
disclose analyzing the sample by matrix-less light desorption/ionization mass spectroscopv, after 
the sample has been applied to the deposited continuous film. Farmer et al. do not disclose 
depositing a continuous film having optical properties and species adsorption properties 
essentially the same as the optical properties and the species adsorption properties of the 
continuous film bulk material and Farmer et al. also do not disclose analyzing the sample by 
matrix-less light desorption/ionization mass spectroscopy, after the sample has been applied to 
the deposited continuous film. Therefore combining Bal et al. v^ith Farmer cannot be used to 
establish nor disclose depositing a continuous film having optical properties and species 
adsorption properties essentially the same as the optical properties and the species adsorption 
properties of the continuous film bulk material and combining Bal et al. with Farmer also cannot 
be used to establish nor disclose analyzing the sample by matrix-less light desorption/ionization 
mass spectroscopy, after the sample has been applied to the deposited continuous film. 

Under a 103 rejection, a prima facie case of obviousness of the invention is made in view 
of the scope and content of the prior art. In order to establish a prima facie case of obviousness, 
"there must be some suggestion or motivation, either in the references themselves or in the 
knowledge generally available to one of ordinary skill in the art, to modify the reference or to 
combine reference teachings. Second, there must be a reasonable expectation of success. 
Finally, the prior art reference (or references) must teach or suggest all of the claim limitations." 
MPEP § 2143. 

Neither Bal nor Farmer separately or in combination teach all the patentable features of 
independent claim 1 and the dependent claims. In light thereof. Applicants respectfully traverse 
the 35 U.S.C. 103 rejection of the claims. 
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In conclusion, in view of the above remarks, Applicants respectfully request the 
Examiner find all claims, 1-3, 6-7, 10-12, 14-19, 66-67 and 69-70 and 125, allowable over the 
prior art and pass this case to issue. 

Since the total number of claims is less than the number of claims already been paid for, 
no additional fees are required. However, if fees are required, they should be charged to Deposit 
Account No. 03-2410, Order No. 30626-101. 

In accordance with Section 714.01 of the MPEP, the following information is presented 
in the event that a call may be deemed desirable by the Examiner: 



ORLANDO LOPEZ (617) 854-4000 



Dated: April 26, 2004 



Respectfully submitted, 

Stephen J. Fonash, et al., Applicants 



By: 




Orlando Lopez / 
Reg. No. 46,880 
Agent for Applicants 
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Background 

Nafion®, was developed by Dr. Walther Grot at DuPont in the late 1960's by modifying Teflon®. Nafion 
was the first synthetic polymer ever developed with ionic properties, and it started an entirely new class 
of polymers called ionomers. These ionic properties of Nafion were created by adding sulfonic acid groups, 
a chemical with very strong ionic properties, into the bulk polymer matrix. Nafion combines the physical 
and chemical properties of its Teflon base material with ionic characteristics that give the final material the 
following properties: 



^ Like Teflon, Nafion is extremely resistant to chemical attack. According to DuPont, 
only metallic alkali metals (sodium in particular) can attack Nafion directly under 
normal conditions of temperature and pressure. This means Nafion does not release 
fragments or degradation products into the surrounding medium. 

\^ Like Teflon, Nafion has relatively high working temperatures compared to many 
polymers. Nafion is used in some applications at temperatures up to 190° C. 

^ Unlike Teflon, Nafion is highly ion-conductive. It functions as a cation exchange 
polymer. 

Nafion is a super-acid catalyst. The sulfonic acid groups attached to the Teflon 
backbone within Nafion function as an extremely strong proton donor due to the 
stabilizing effect of the large polymer matrix attached to the sulfonic acid. 

Q Nafion is very selectively and highly permeable to water. The sulfonic acid groups in 
Nafion have a very high water-of-hydration, so they very efficiently absorb water. 
Interconnections between the sulfonic acid groups lead to very rapid transfer of water 
through the Nafion. 



Applications 

These unusual properties are exploited in the four main applications of Nafion today. 



I Nafion is used to fabricate ion-exchange membranes used to produce chlorine gas 
and lye (sodium hydroxide) by the electrolysis of salt water. 

2. Nafion is used to selectively dry or to humidify gases. The largest segment of this 
application involves drying or humidifying breath for anesthesia, respiratory care, or 
plethysmography, but many applications exist in the laboratory and in industrial use. 



3^ Nafion is used as the proton exchange membranes in polymer electrode fuel cells. 



http://www.perniapure.com/ncwweb/Nafion%20propcrties.htm 
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4 Nafion is used as a super-acid catalyst in the production of fine chemicals. ; 
Physical Appearance 

■I Nafion is initially produced in a salt form with its sulfonic acid groups neutralized. In this form, Nafion is I 
thermoplastic (extrudable and formable using heat and pressure), but not chemically active. Once j 
extruded into its final form (typically sheet or tubing), Nafion is chemically activated by conversion of the 
salt form to the acid form. At this stage, Nafion is a translucent plastic that is similar in appearance to 
Teflon, but clearer and less opalescent. 

Once activated, Nafion immediately begins to react with its surrounding environment. Moisture is ' 
absorbed and exchanged with the surroundings. Reactive organic gas components within the surrounding i 
air may undergo chemical changes when exposed to the super-acid catalytic activity of Nafion. These i 
; organic gases may combine to form larger compounds that are liquid or solid in nature. Over time, these j 
; organic residues build up a deposit on the Nafion. The original color of the Nafion changes gradually from , 
I translucent to yellow, then brown, then even black. ' 

!i 

ii The chemical nature of the Nafion itself is not changed by this process. The Nafion is not directly 

attacked, and the Nafion is not degraded. Its properties of water permeability, ion-exchange, and acid 
il catalysis are not directly affected. l| 

For Nafion to transport water (dry or humidify gases), the water must be able to reach the active sites on 

the surface of the Nafion. In extreme cases, deposits of organic residues can build up on the Nafion to 
il the point where the surface of the Nafion is sufficiently occluded to reduce its functionality. When this | 

occurs, the discoloration of the Nafion is severe, black rather than merely yellow or brown. Even in these 
I extreme cases, Nafion is still functional, it merely loses efficiency. Loss of a few percent of its total ] 
^ functionality due to partial occlusion means that a gradual decline in performance will be observed, not a ] 

catastrophic failure. 

Under typical storage or operating conditions, Nafion will turn somewhat yellow within a year and brown i 
within three to five years. Even after turning brown, the Nafion is still fully functional, and no loss of i 
I performance efficiency is likely to be observed. | 

These unwanted chemical reactions in the air surrounding Nafion are stimulated by exposure to light and i 
[ to elevated temperatures. Storage in sealed bags in the dark will extend the pristine original appearance | 
of Nafion. J 

The original clear appearance of Nafion can be restored by cleaning. Various cleaning solvents may be ij 

I used if desired, usually non-polar solvents such as hexane. The most effective cleaning method (used I 
:1 by Perma Pure Inc. in its dryer service procedures) is to boil it in strong acid. Nafion is not damaged by 

;i these procedures, but the discolored surface residues are eliminated. i 
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Crystallinity in Cast Nafion 

M. Ludvigsson, J. Lindgren/ and J. Tegenfeldt* 

The Angsirrym Lcboratoiy, Inoj-ganic Chcmisfiy. Vppsah Unhersity. 5E-75! 21 Uppsala. Sweden 



The pcrfomiancc of a polymer clcch olyfc fuel cell is cntically dependent on ihc watcf uptake in the polyincr electrolyte, usually 
Nafion. Nafion in solmion isi often painted onto Ihc cIcc»Todcs of Ihc fuel cell. It is iinportani thai this cast Nafion film' stay amor- 
phous and not crvsiollizc. Cast Nafion films, ca. 1 |j.ni thick, crystallized on silicon plates when kept in air ai room temperature Tor 
a long time. The film.s contain large ciyslallinc regions ranging from 0.5 mm to several millimeters in size. X-ray diffraction (XRD) 
and Fourier rransfomi infrared (FTIR) microspcciroscopy have been used to investigate the ciysialline and amoi-phou.s regions. The 
XRD shows TWO sharp peaks. One of the pcnks is dcx clopcd before die second one appears in ihc diffracfogram, indicating that 
there might be Iwo lypcs of crystallizing processes. FTIR spectra of ihc amoi-phous and crysiallinc regions in ihc films show impor- 
tant differences. In ihc cry^stallino regions, the film contains the sulfonic acid at the end of the .side chains: hence, the crysiallinc 
rcgionjs contain no water molecules. In the amorphous regions there is a complete proton transfer from the acid to the water mol- 
ecules, and sulfonate groups arc obtained. 

2000 The Elccirochcmical Society. SOO 1 3-465 1(99106-021 -8. All rights reserved- 



Manuscript submillcd .lunc 2. 1999; revised manuscript received December 10, 1990. 



Nafion^^^ is a pcrfluoTosuIfonic polymer membrane having ex- 
ceptional chemical, thermal, and mechanical stability. It can act as a 
cation exchange membrane and it has found application as a separa- 
tor in chlor-alkali cells. More recently, and due to its high proton 
conductivity, Nafion has attracted much interest as a polymer elec- 
trolyte membrane in polymer electrolyte fuel cells (PEFCs). Rcccm 
reviews summarize the applications as well as the structure and 
properties of perftuorinated membranes. ^--^ To enhance the perfor- 
mance of the PEFC, Nafion in solmion is painted onto the electrode 
surfaces. Results from wide-angle X-ray diffi-action (WAXS) studies 
arc rather inconsistent and imply That there is either no crystallinity 
in Nafion membranes,-' that they are partly crystalline.^ or that cry.s- 
tallinity- can be introduced in recast amorphous films by thermal 
treatment.^ An electron microscope study suggests that crystalline 
Nafion has a linear zigzag structure Just as polyethylene.^ This poly- 
ethylene-like eirucrure is also supported by an investigation of ther- 
mally treated recast Nafion.' 

Our question was whether it would be possible to introduce crys- 
tallinitv' into Nafion with some sort of hca\'y ion irradiation. We de- 
cided on xenon ions. Cast Nafion films iiTadiaced in vacuum by xenon 
ions were in fact observed to crystallize. Even though the films were 
exposed to different doses of ions, they crystallized in the same way. 
The size of the crystalline regions was 0.5-10 mm and thus easily 
observed. This attracted our attention and some further investigations 
on this material were made, as well as on films prepared without any 
irradiation. 

Experimental 

Thin films (ca. I |xm thick) of Nafion (cquiv wt I 100) were cast 
on single-crystal silicon windows from a solution of lower aliphatic 
alcohols and water, Fluka Chcmika. For use in another experiment, in 
an ongoing project, the films were put in vacuum, I0~^ mbar, for 
approximately 1 0 h while they were irradiated with xenon ions. The 
irradiation was made on the cyclotron facility at the The Svcdberg 
Laboratoiy, Uppsala University. Each film was exposed to between 1 
and 100 ions-'p.m^. The films were siored at room lemperamre in the 
laboratory for future investigations. 

After some time, large crystalline regions. 0.5-10 mm, were fomi- 
ed in the films (Fig. 1 ). If a drop of water was added to the crystalline 
region of the Nafion film, it became amorphous immediately. Fouri- 
er transform infrared (FTIR) microspectroscopy (reflection mode) 
and X-ray diffraction (XRD) were used to investigate the aniorphous 
and ciy stailine regions. Due to the size of the crystalline regions, the 

* Ekcirochcinical Society Aclivc Member, 
' E-nnuil: .1iin.LindiJiPcn(£i;tcemi.uu.so 



infrared beam could be focused on cither crystalline or amorphous 
regions to record infrared spectra in open oir. 

Funhcrmore, in another ongoing project, reference spectra at dif- 
ferent relative humidities (RH) were recorded.'^ This experiment 
included a cell where the desired RH and temperamrc could be 
achieved. In the experimental arrangement it was possible to achieve 
a state where n = \ (n - H2O/SO3 ). i.e., in \^cuum at room tem- 
perature,^ The ;i = 0 state was achieved in vacuum at an elevated 
temperature, .'\ Nafion film saturated with n = 14 could be obtained 
in i00%RI-I.'<> 

The infrared spectra were recorded on a Bio-Rad/Digilab FTS 
45 FTIR spectrometer (equipped with a microscope) covering 400- 
4400 cm" ' with a resolution of 2 cm"'. The X-ray diffracto grams 
were recorded on a Siemens D5000 diffractomctcr with Cu Ka 
radiation. 

Results and Discussion 

In Fig. 2, the infrared spectra of the crystalline and amorphous 
regions of a xenon-irradiated Nafion film are shown together with 
spectra of Nafion films containing, 0. 1 , and 14 water molecules per 
sulfonate group.'^ The symmetric stretching band of the SO3 -group • 
tenninating the side chains of Nafion is situated at ca, I06O cm~L 
In the top spectrum (a) of crystalline Nafion two bands at ca. 1400 
and 910 cm"' appear. These two bands correspond to the v(S=0) 
and v(S-OI-l) .stretching modes of the SO;,H group. These two bands 
do not appear in the IR spectrum (b) of the amorphous region of an 




Figure I, CtyJi tall ires in Nafion exposed to high vacuum and xenon ions. 



Received Apr-05-04 OSillpm From- 



To-PERKINS SMITH Page 02 



Rpr- 05 04 02: 12p 



1304 



Journal of The Elearochemicu} Socier\\ 147 (4) 1303-1305 (2000) 
S00I3-465 1(99)06-021-8 CCC: S7.0O i.'i Tlic rjcctrochemical Society, Inc. 




1400 



1200 



1000 



wavenumbers (cm' ) 



Figure 2. IR spcctiti of (j1 the co'stallinc region and (b) the amorphous 
region of a ci>iitalline Nafion film exposed to high vacuum and xenon ions 
and Nnfion films contiiininii (c) 0. (d) 1, and (c) 14 water molecules pcv sul- 
fonate group. 



irradiated film. They arc present in a completely dried film as in 
spectrum (c) where there is no water at all. However, as soon as there 
is more Than one water molecule per sulfonate group, the hydrogen 
is transferred to the water molecule, i.e.. the two bands arc missing 
in spectrum (d) and (e).*^ This indicates that the side chains are ter- 
minated by SO3H groups in the dry crystalline region and not by 
SO;,' groups as in the amorphous region. Thus, the cr^'stallinc rc; 
gions contain less water than the amorphous regions. This could be 
an additional reason why Nafion membranes absorb less water at 
elevated Temperatures, provided that rhe Nafion membrane crystal- 
lizes partly at elevated temperatures." 

[n the diffractograms (Fig. 3>. the crystalline peaks of the Nafion 
film appear on top of the amorphous band. There arc two broad amor- 
phous features corresponding to d-spacings of 5.0 and 2.2 A. In The 
upper diffractogram of a xenon irradiated Nafion there arc four cr\'3- 
tallinc peaks, see Table T. Tn the lower diffractogram of a Nafion film 
kept at room temperature in the laboratoi-y. only one peak is obscn'cd. 

Porat et ui have made an electron microscope investigation of 
very thin Nation films, 30-60 nm (compared to ca. 1 p-m in the pre- 
sent study)-*" They found similarities with the electron microscopy 
data of polyethylene, sec Table II. They concluded thai the crys- 
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Figure 3. Diffi-actograms of I ^im thick Nofion films; (top) filin containing 
easily ohscrvabic cr\'staninc regions (exposed to high vacuum and xenon ion 
bombardment) and (bonom) film kept ni room tcmpcraiurc and normal 
ntmo spheric conditions for A months. 



Table I. The d-spacings of (he cr\'stallinc Nafion film in Fig. J. 
Calculated values arc obtained from a least-squares fit. 
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Intensity 
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4.21 
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3.32 


3.32 
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2.72 
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Strong 
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2.10 


2.10 


Weak" 



Table II. The d-spacin«s and suggested lattice parameters for 
cr>'stalltnc polyethylene and >afion.*-'* 



d-Kp^cing 
(A) PE 



d-spticing (A) 
Porat L>i al. 



fllO) 
(200) 
(020) 
(310) 



4.1 
3.7 
2.4 
2.2 



A.]5 
3.77 
2.A9 
2.24 



talline regions from the tluorocarbon backbone in their Nafion films 
had a linear zigzag structure just as polyethylene, '"^ and not a twist- 
ed chain as polytctrafluoroethylene (PTFE).'^ The crystalline re- 
gions in the present study have sphenilitic morphology. In general, 
spherulites grown from a bulk sample consist of lamellar cr>'stalline 
plates of folded polyrncr chains which arc perpendicular to the plane 
of growth, i.e., the lamellae and the film itself (Fig. 4).'**-'^ This is 
equivalent to a c axis oriented perpendicular lo the plane of the 
Nafion film, and it can thus be assumed that we have the same pre- 
ferred c axis orientation normal to the plane of the film as in the 
study by Porat er at. Thus, it is possible to calculate the lattice para- 
meters a and h ustna the expression for an onhorhombic system: 
d"- ^ {hla)^ 4- ikih)'^ + WcY^ The result is a = 6.64 ± 0.0 lA and 
^> = 5.44 r 0.01 A (for Porat et al, a = 7.49 A and b = 5.02 A). It 
is not possible to determine the lattice parameter c since part of the 
necessary* information is missing. Our results thus indicate an or- 
thorhombic cell but with different cell parameters from those of 
Porat ej al. for Nafion and for polyethylene. 

A sample of a cast Nafion film was kept in the laboratory' at room » 
temperature and under noiTnal atmospheric condition.*; for an extend- 
ed period of time to see if the crj'siallization is a spontaneous process 
or if the treatment with high vacuum and xenon ions is the cause. .\ 
diffractogram recorded directly after the film was cast shows only the 
broad amorphous features in Fig. 3. After a week, however, the strong 
peak at d = 2.72 A appears in the diffractogram (Fig. 3). but 4 months 
larer nothing else had happened and there was still only one er>*s- 
talline peak. This may indicate that the full crystallization does not 
occur unless initiated, e.g.. by irradiation of xenon ions or by appli- 




Figui-o 4, Schematic drawing of lamellar plates in the ciystallinc regions of 
sphci-ulitic morphology. 
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cation of high vacuum and that some scmico'stallinc phase develops 
in a slow process in this film without the initiation. The exact coinci- 
dence of the single observed diffraction peak from the nonirradiated 
sample with the (020) peak from the irradiated sample would imply 
an ordering in one dimension in the b direction of the same type in 
both samples. The fact that the ci-ystalline region became amorphous 
when wetted raises the qxiesiion whether the irradiated poKnner is 
crystalline because it is dried out or dehydrated because it is crys- 
taUine? Since only some regions of rhe Nafion films are crystalline 
and large amorphous areas are observed, we think that the latter state- 
ment is tmc. Otlicrwisc the films ought to have crystallized com- 
pletely. The introduction of these crystalline regions will probably not 
be of any harm in a fiiel-eell application because they have been forc- 
ed there with the help of xenon ion bombardmcnr in high \^cuum. 
Nevertheless, if they would appear in a dried-out membrane in an 
operating fuel cell, they arc reverted to amorphous regions as soon as 
it is hydrated again. The silicon substrates seem to have no direct 
effect on the crystalline regions, since no preferred orientation can be 
found on any of the substrates. 

Conclusions 

In this work we used XRD and IK spectroscopy to study the crys- 
talline phase of cast Nafion filnis. The earlier reported polyethylene- 
like orthorhomhic structure of thin Nafion films is not supported by 
our results that indicate an orthorhombic strucnire with different cell 
parameters. There is, unfortunately, not enough information in the 
results from the XRD for a conclusion regarding the exact chain con- 
fomiation or packing at present. Ncvcnhclcss, the IR spectra imply 
that the crjr'stalUnc regions of the Kafion films contain less water 
than the amorphous regions. 
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Nation® is a commercially available perfluorosulphonate cation exchange membrane 
commonly used as a perm-selective separator in chlor-alkali electrolysers and as the 
electrolyte in solid polymer fuel cells. This usage arises because of its high mechanical, 
thermal and chemical stability coupled with its high conductivity and ionic selectivity, 
which depend strongly on the water content. The membrane was therefore studied in 
different states of hydration with two complementary techniques: atomic force microscopy 
(AFM) and small angle X-ray scattering (SAXS) combined with a maximum entropy 
(MaxEnt) reconstruction. Tapping mode phase imaging was successfully used to identify 
the hydrophobic and hydrophilic regions of Nafion. The images support the MaxEnt 
interpretation of a cluster model of ionic aggregation, with spacings between individual 
clusters ranging from 3 to 5 nm, aggregating to form cluster agglomerates with sizes from 
5 to 30 nm. Both techniques indicate that the number density of ionic clusters changes as a 
function of water content, and this explains why the bulk volumetric swelling in water is 
observed to be significantly less than the swelling inferred from scattering measurements. 
® 2000 Kluwer Academic Publishers 



1. Introduction 

Nafion® is a commercially available perfluoro- 
sulphonate cation exchange membrane (CEM) manu- 
factured by E I du Pont de Nemours & Co. Inc. It is 
commonly used as a permselective separator in chlor- 
alkali electrolysers [1,2] and as the electrolyte in solid 
polymer fuel cells (SPEC). Perfluorosulphonate cation 
exchange membranes are used in these applications be- 
cause of their high ionic conductivity and their high me- 
chanical, thermal and chemical stability. Structurally, 
Nafion consists of a hydrophobic tetrafluoroethylene 
(TFE) backbone with pendant side chains of perfluori- 
nated vinyl ethers terminated by ion-exchange groups. 
The chemical structure of Nafion is shown in Fig. 1. 
The ion content can be varied by changing the ratio of 
the two components. 

Other perfluorosulphonate cation exchange mem- 
branes with similar structures have also been developed 
by the Asahi Chemical Company (Aciplex®) and the 
Asahi Glass Company (Flemion®). The Dow Chemical 
Company also developed a material with a shorter side- 



chain than those of Nafion and the other perfluoro- 
sulphonaies [3). 

The ionic conductivity in perfluorosulphonate mcnv 
branes is important since ii should be as high as possible 
to minimize ohmic losses in chlor-alkali electrolysers 
and to maintain high output power densities in fuel cell 
applications. The membrane conductivity is strongly 
influenced by the water content [4]. Water manage- 
ment is also an important problem for efficient SPFC 
operation to avoid flooding of the gas diffusion elec- 
trodes. A number of factors affect the water content 
including the cation form, ion exchange capacity of 
the membrane and equivalent weight (EW). The ion 
content is usually expressed in terms of the equivalent 
weight of the polymer. The equivalent weight is defined 
as the weight of dry polymer in grams containing one 
mole of exchange sites. The desired equivalent weight 
is achieved by varying the ratio of vinyl ether monomer 
to TFE. The useful range of equivalent weights is from 
600-1500 g per mole, between the solubility and perco- 
lation limits. The industrial applications of Nafion have 
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prompted considerable research effort summarised by 
Eisenberg and Yeager in 1982 [5] and more recently by 
Tant etaL in 1997 [3]. 

Perfluorosulphonate polymers typically have ordered 
structures in which the hydrophilic end groups aggre- 
gate within a hydrophobic matrix composed of the fluo- 
rocarbon backbone of the polymer. The density contrast 
between the ionic aggregates and the matrix gives rise to 
scattering of X-rays [6-9] and neutrons [5]. The diffrac- 
tion patterns show a broad ring reflection, the so-called 
'cluster' reflection, with a peak position of 3.5-5.5 nm, 
together with a scattering upturn at low angles. Many 
different structural models have been proposed to ex- 
plain these features, and they can be divided into two 
main groups. The first type are the interparticle models, 
such as the two-phase model [5, 10] and the lamellar 
model [5, 1 1 ], in which the scattering is produced by the 
interference between ionic aggregates. The second are 
the intraparticle models, such as the core-shell model 
[5, 12, 13], in which the scattering is due to fine struc- 
ture within an individual ionic aggregate. 

There have been a number of electron microscopy 
studies of Nafion [14, 15] which support the clus- 
ter model of phase separation, although the size scale 
of the structures observed did not correspond directly 
with those detected by scattering studies. The Cluster- 
Network model [5, 16] postulates a large scale organi- 
sation of clusters with transient connective tubes which 
are in constant flux, A fully-reversible dynamic reor- 
ganisation of the clusters is proposed to occur on re- 
hydration. Positron annihilation spectroscopy has also 
been used to study different cation forms of Nafion 
membranes by measuring the free volume. These stud- 
ies [17-19] found that the formation and expansion 
of clusters is always associated with a change in 
free-volume structure resulting in smaller free-volume 
holes. 

There have been two recent studies undertaken using 
scanning probe microscopy (SPM) [20, 21], to investi- 
gate the structure of Nafion in differing states of hydra- 
tion. A number of swelling studies have been carried out 
on Nafion [22, 23] using a wide variety of solvents with 
some remarkable results, including a swelling of 360% 
with tributylphosphate (TBP) as used by Lehmani et ai 
in their SPM study [21]. Nafion has proved to be a suit- 
able sample for other forms of SPM including scanning 
tunnelling microscopy (STM) and scanning electro- 



chemical microscopy (SECM) [24, 25]. Fan and Bard 
[24] studied a spin-coated Nafion film using SECM and 
found a domain-like structure, consisting of circular 
structures, 1-2 nm in diameter made up from a con- 
ductive central zone surrounded by a less conductive 
region. 

Since many of the useful properties of Nafion mem- 
branes depend strongly on the water content, studies in 
different states of hydration using the two complemen- 
tary techniques of atomic force microscopy (AFM) and 
small angle X-ray scattering (SAXS), might be infor- 
mative. The difference in probe-specimen adhesion be- 
tween the hydrophobic backbone and hydrophilic side 
group regions of Nafion offers the possibility of ob- 
serving the spatial distribution of these two regions at 
the surface of the polymer using lapping mode phase 
imaging. Assuming that these are representative of the 
bulk, then the number density and average cluster size 
can be deduced as a function of water content, and the 
results compared with those from SAXS experiments. 

2. Experimental method 

2.1. Variation of the mass of Nafion 
with humidity 

The uptake of water in Nafion membranes has been 
studied as a function of environmental humidity using a 
dynamic vapour sorption study (DVS) [26]. The sample 
was suspended in a 0% RH environment and weighed. 
The humidity was increased in steps using different sat- 
urated salt solutions, the sample was allowed to equili- 
brate and then weighed, the process was then repeated. 
A graph of percentage mass gain versus humidity was 
plotted from the data. Given that the equivalent weight 
of the sample and the theoretical weight are known, 
the average number of water molecules per sulphonic 
acid group can be calculated at 0% RH, by comparing 
the theoretical weight to the measured value. If it is 
assumed that all of the water is incorporated into the 
hydrophilic regions, a graph of average number of wa- 
ter molecules per sulphonic acid group can be plotted 
over the full range of humidities. 

2.2. Sample preparation 

All of the experiments were carried out using the most 
commonly used form of Nafion, 1I5H"^, which has 
the following physical properties; a nominal equiva- 
lent weight of 1100 (actually 1070) and thickness of 
5 thousandths of an inch (^127 ^m). 

The membrane will readily rehydrate if left exposed 
to a high relative humidity environment [27]. Careless 
handling can result in the membrane ion-exchanging 
from the acid (H"*") to a salt form (e.g. Na"^ or K+). All 
samples were therefore routinely prepared by refluxing 
with a 50/50 mixture by volume of concentrated nitric 
acid and de-ionized water, then de-ionized water alone, 
to ensure that the membrane is in the H"*" form and free 
from any chemical impurities. 

2.3. Small-angle X-ray scattering 

Two dimensional, point collimated SAXS data were 
collected using nickel-filtered Cu K„ radiation on a flat 
plate Rigaku-Denki camera, with a typical sample to 
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film distance of around 250 mm. The X-ray genera- 
tor was an Elliott GX21 rotating anode generator The 
films were scanned in transmission using an Optronics 
P2000 drum densitometer. The resulting images were 
converted to 8-bit binary files, with a pixel dimension 
of 100 X 100 ^m. The low-angle limit due to beam di- 
vergence was approximately 0.026 nm"*, which cor- 
responds to a maximum discemible size for features 
in real space of around 40 nm. The diffraction pat- 
terns were corrected for Lorentz-polarisation effects 
and sample absorption. 

Previous studies of SAXS from hydrated ionomer 
membranes used heat-sealed polymer bags to enclose 
the samples during exposure [6, 7]. This method was 
found to be unsatisfactory for quantitative work due to 
fluctuations in the environmental humidity level during 
the long periods (up to 24 hours) required to collect 
the SAXS data. An alternative approach, based on the 
controlled humidification of air using concentrated salt 
solutions, was adopted [28]. 

2.4. Maximum entropy method 

A full account of the maximum entropy (MaxEnt) 
method used to interpret SAXS data has already been 
published [29], and only a summary will be presented 
here. The aim of the MaxEnt technique is to reconstruct 
a two dimensional charge density map which is a good 
fit to the observed scattering data. The MaxEnt criterion 
is required in order to resolve the inherent ambiguity 
in this fitting process due to the absence of phase infor- 
mation in the experimental data. Since there are many 
charge distributions which can give rise to identical 
diffraction patterns, we select the one with the highest 
entropy for a given goodness-of-fit. Of course, there is 
no guarantee that this model will necessarily be stereo- 
chemically or thermodynamically realistic, but if a rea- 
sonable structural interpretation can be found then we 
can be sure that this it is the most statistically likely in- 
terpretation which is consistent with the scattering data. 
In this way, the MaxEnt method is superior to other ap- 
proaches which assume a priori a structural model. The 
algorithm used to carry out the reconstruction process 
is based on the ^Cambridge Algorithm' [30]. 

2.5. Atomic force microscopy 

There are two principal modes of atomic force mi- 
croscopy (AFM) [31] operation: contact mode and tap- 
ping mode. In contact mode, the tip remains in close 
contact with the sample, operating within the repulsive 
regime. In tapping mode, the cantilever is oscillated 
near to its resonant frequency, above the surface, only 
contacting the surface briefly during each cycle of its 
oscillation. Tapping mode is therefore more suitable 
for imaging delicate samples owing to the lower lat- 
eral forces and has been applied to many polymer sys- 
tems [32-35]. It also has the added advantage of being 
able to obtain phase images as well as topographical 
data. Tapping mode phase imaging is a relatively new 
AFM technique, and it can differentiate between areas 
with different properties regardless of their topograph- 
ical nature [36-42], The phase angle is defined as the 
phase lag of the cantilever oscillation relative to the 



signal sent to the piezo driving the caniilever. Us value 
depends on the energy dissipated in the lapping inter- 
action of probe and specimen [43, 44]. 

Circular samples 10 mm in diameter were mounted 
on magnetic stainless steel stubs, placed in the Digital 
Instruments Extended Multimode AFM and imaged in 
tapping mode with silicon cantilevers in order to pro- 
vide topographic and corresponding phase images. The 
membrane and microscope were placed in a purpose- 
built environmental chamber to control the humidity 
and the samples were dehydrated overnight. A variety 
of techniques was employed to reduce the humidity 
including silica gel, phosphorus pentoxide, a liquid ni- 
trogen cold finger and dry nitrogen gas. A humidity 
minimum of (0.8 ± 2.0)% was obtained using nitrogen 
gas passed through molecular sieve material. The sam- 
ple was then imaged and allowed to rehydrate slowly to 
room humidity (34 ± 2)% over the course of 12 hours 
whilst being continuously imaged, and the humidity 
recorded. 

2.6. Cluster counting algorithm 

Phase images from the beginning and end of the rehy- 
dration sequence were analysed with the aid of a clus- 
ter counting algorithm [28, 45]. The 256 x 256 images 
were exported, cropped, saved as eight-bit data and then 
converted to binary. A series of thresholds from 0 to 
256 was then applied to all the data points in the image; 
if the point was equal to or greater than the threshold 
value it was set to one (white), and, if not, it was set to 
zero (black). A cluster is defined as an isolated group 
of white pixels. The number of unique clusters, total 
area of clusters and average cluster area was then cal- 
culated for each threshold. The number of clusters at 
each threshold value was plotted for the two humidity 
extremes. The same algorithm was also used to analyse 
the MaxEnt reconstructions of the SAXS data acquired 
at ambient and 100% relative humidity. 

3. Results and discussion 
3.1. Variation of the mass of Nation 
with humidity 

Analysis of the DVS results is more revealing than pre- 
vious studies [27]. In addition to a graph of percent- 
age mass change against humidity ( Fig. 2i), a graph of 
the ratio of water molecules to sul phonic acid groups 
against humidity (and not just the additional number) 
can be plotted (Fig. 2ii). This has been calculated us- 
ing a theoretical dry mass for Nafion and by assuming 
that the EW is 1070. The difference in mass between the 
theoretical dry mass and that at 0% RH is the amount of 
water bound to the sulphonic acid groups. The average 
number of water molecules per sulphonic acid group at 
0% RH was found to be 1.55. In the same regime as 
that used for the AFM data, in the humidity range of 
(0-50) ± 2%, on average an additional water molecule 
was added to each sulphonic acid group for an increase 
in RH of 14 ±2%. At higher humidities this becomes 
much more rapid and the assumption that all of the wa- 
ter is incorporated into the hydrophilic regions become 
less valid as water will start to condense on all regions 
of the membrane. 
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Figure 2 Percentage mass gain (i) and Number of addition water mole- 
cules per sulphonic acid group (ii) against relative humidity for Nafion 
115 H"^. The relative humidity was cycled from to 0 to 100% and back 
again. 
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Figure 3 Peak position of the 'cluster' reflection from Nafion 1 15 H*^ 
membrane as a function of the volimie fraction of absorbed water. 

3.2. Small-angle X-ray scattering 

The variation in the peak position of the ^cluster* reflec- 
tion in the SAXS from a Nafion 1 15 H"*" membrane as 
a function of the volume fraction of absorbed water is 
shown in Fig. 3, along with the data published by Gierke 
et al. [6] for comparison. There are no values below a 
volume fraction of 0. 1 as it was not possible to detect 
the 'cluster' reflection at such low water contents due 
to lack of electron density contrast. The dependence of 
the peak position on water content is predicted by the 
type of scattering model which is used to interpret the 
data. If the ^cluster* reflection were due to the coher- 
ence of intercluster spacings [5, 10], or the incoherent 
sum of interlamellar spacings [5, 11] then it might be 
expected that the increase in peak position would be 
directly proportional to the volume of absorbed water. 
However, if the reflection were produced by the scatter- 
ing from individual three dimensional clusters models 
[5, 12, 13] then its spacing should vary with the cube 
root of the water volume. In fact, neither of these two 
simple behaviours occurs. 



In order lo explain the observed swelling behaviour, 
it is instructive lo consider a MaxEni inierpreiaiion of 
the scattering data. Fig. 4 shows MaxEni reconsiruc- 
tions from membranes under ambient and 100% RH 
conditions. The MaxEnt reconstructions should be in- 
terpreted as 2-dimensional projections of the electron 
density distributions within a representative volume of 
membrane. Regions of low electron density appear light 
grey or white, while regions of high electron density ap- 
pear dark grey or black. The absolute electron density 
range in each case is scaled to fit the grey-scale. The re- 
constructions are not unique, belonging lo a large set of 
possible maximum entropy solutions, all of which arc 
capable of reproducing the experimental data subject to 
the usual errors associated with X-ray counting statis- 
tics. However, the region of membrane represented in 
the reconstructions is sufficiently large that it is able to 
reproduce a wide range of structural variation within a 
single electron density map. 

Both MaxEnt reconstructions contain fragmented re- 
gions of low charge density with a broad range of sizes. 
There is an increase in electron density contrast be- 
tween these regions and the background, along with a 
coarsening of the size distribution from Fig. 4i (ambi- 
ent) to Fig.4ii (100% RH). In order to make the physical 
interpretation of these features clear, the reconstruc- 
tions were spatially filtered to separate the small scale 
structures associated with the 'cluster' reflection from 
the large scale structures associated with the scattering 
upturn. The filtering was performed by windowing out 
the unwanted regions of the diffraction patterns, prior 
to MaxEnt reconstruction. The high spatial frequency 
images, shown in Fig. 4iii and iv for ambient and 1 00% 
RH respectively, both contain collections of point-like 
scatterers. However, as the water content is increased 
from Fig. 4iii to Fig 4iv, the spacing between the scat- 
terers increases, and their number density decreases. It 
is the coherence of the spacings between the scatter- 
ers that is responsible for the 'cluster' reflection, and 
this interpretation is borne out by reconstructions from 
highly oriented samples [45], We therefore identify the 
point-Uke scatterers as being individual ionic clusters. 

The low spatial frequency images, shown in Fig. 4v 
and vi, show large scale structures which are composed 
of agglomerates of the clusters identified in the high fre- 
quency images. Although the agglomerates do not swell 
appreciably with increasing water content, the electron 
density contrast between these features and the poly- 
mer matrix increases noticeably due to the difference 
in density between the water and fluorocarbon phases. 
The shape of these cluster agglomerates, which is re- 
sponsible for the form of the low angle upturn in the 
SAXS pattern, is related to the spatial coherence of the 
clusters. We therefore predict an intimate relation be- 
tween the form of the uptum and the degree of arcing of 
the *cluster' reflection. If the MaxEnt interpretation of 
the SAXS data is correct, then it should not be possible 
to see an arced 'cluster' reflection unless the low angle 
uptum is also anisotropic. 

Assuming the MaxEnt interpretation is correct, then 
it is useful to calculate the relative number density of 
ionic clusters in Fig. 4iii and iv and relate this to the dis- 
crepancy between the bulk swelling of the membrane 
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Figure 4 MaxEnt reconstruction of SAXS data from Nafion 115 H+ membrane: (i) ambient humidity and (ii) 100% RH. The reconstructions have 
been high pass filtered in (iii) and (iv) and low pass filtered in (v) and (vi). 



and the swelling which is inferred from the increase in 
the peak position of the *cluster* reflection. Although 
the mean lateral intercluster separation increased by 
44.1%, the overall lateral expansion of the membrane 
is expected to be significantly less than this due to 
a decrease in the number density of clusters. Unfor- 
tunately, it is difficult to predict the resultant macro- 
scopic swelling quantitatively, because the electron 
density distributions are a 2-dimensional representation 
of 3-dimensional space. It is also possible that system- 
atic changes in the spatial ordering of the clusters on 
swelling need to be incorporated into the model. Nev- 
ertheless, the mechanism described is qualitatively ap- 
pealing and has been derived with the use of only very 
general a priori hypotheses, i.e.: (i) the use of a Max- 



Ent electron density distribution, (ii) the legitimacy of 
spatial filtering and (iii) the definition of a 'cluster' as 
an isolated agglomerate of pixels. 

Although the redistribution of exchange groups was 
first predicted by the Marx infinite paracrysialline 
model [10] it is significant thai this phenomenon can 
be inferred directly from the small angle X-ray data. 
In addition, there is also evidence for such a redistri- 
bution from phase contrast atomic force microscopy as 
discussed in the following section. 

3.3. Atomic force microscopy 

Six 1 ^m phase images from a sequence of 56 images 
of 1070EW Nafion obtained over a range of humidities 
from (9-34) ± 2% are shown in Fig. 5i-vi. The phase 
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Figure 5 Six 1 fim tapping mode phase images of the same area of Nafion at relative humidities of (i) 9, (ii) 13, (iii) 19, (iv) 23, (v) 28 & (vi) 33 ± 
2% (Z scale = 30°). 
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contrast between the hydrophilic sulphonic acid groups 
and the hydrophobic regions increases with humidity 
due to preferential sorption of water to the hydrophilic 
regions. It is very difficult to monitor the same area 
as the membrane expands both vertically and latera lly. 
The same feature has been circled in each of the images 
to aid comparison and to illustrate the drift to the right 
which occurred throughout the series. 

Cluster-like structures with a diameter of 5-30 nm 
are clearly visible, particularly in the phase images. 
These structures do not always correspond to features 
in the height image. Features at the same height can 
have a different phase signal and vice versa, indicating 
minimal topographic coupling. This is consistent with 
the earlier AFM studies [20, 21]. Chomakova-Haefke 
et al. [20], produced samples by evaporating a dilute so- 
lution of Nation on to gold-sputtered cover slips. They 
were found to have a fibril network structure, which be- 
came more ordered as the membrane swelled. Lehmani 
et al. [21], imaged commercially available nominally 
1 100 EW Nafion 117. It was dehydrated in a vacuum 
oven at 80 °C and was first imaged "dry", then swollen 
with de-ionized water, and finally swollen with TBP. A 
super-structure of spherical domains with an average 
diameter of 45 nm containing 1 1 nm grains was found. 
According to the section analysis of the microstructure, 
the interstitial regions have a mean thickness of 5 nm, 
which may correspond to the cluster size. 

The corresponding topography images from the be- 
ginning and end of the rehydration sequence are shown 
in Fig. 6i and ii respectively. The same area has been 
circled in both images demonstrating the ability to im- 
age the same area whilst following a dynamic process. 
There is no appreciable change between the two im- 
ages, therefore any topographic coupling in the phase 
images remains constant and is not responsible for the 
increase in phase contrast. 

Cluster counting analysis of the phase images from 
the beginning and end of the rehydration sequence is 
shown in Fig. 7i and ii, respectively, and tabulated in 
Table I. The sharp peak in Fig. 7i implies that all of 



TABLE I The number of clusiers. touil percentage area of cluster and 
average cluster size for I /xm wide phase images of Nahon al (9 & 34) 
±2% relative humidity 

Humidity ±2% 9 34 

Number of Cluster 1415 1 1 18 

Total Percentage Cluster Area 1 8.9 1 7.6 

Average Cluster Area (nm-) 1 34 1 58 




0 50 100 150 200 250 300 

Threshold 




0 50 100 150 200 250 300 



Threshold 

Figure 7 A graph of number of clusters against threshold value for the 
two humidity extremes of 9 & 34 ± 2%. 

the clusters have a similar phase angle and iherefore 
energy loss at lower humidities. At higher hutnidities 
as shown in Fig. 7ii, the peak has moved lo higher 
thresholds as the average energy loss has increased; 
the peak has also broadened significantly as the clus- 
ters have a wider range of energy losses. The average 
cluster area of between 134 and 158 nm^ is consistent 
with the images shown in Fig. 5, where clusters have 
a range of linear dimensions between 5 and 30 nm (or 





Figure 6 Topography images corresponding to those in Rg. 5i and vi for the humidity extrenrres of 9 & 34 ±2%. (Z scale = 50 nm). 
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approximate areas of 25 to 900 nm^). This is signif- 
icantly larger than the individual clusters which give 
rise to the ^cluster* reflection, but of comparable size to 
the cluster agglomerates. 

The decreasing number of clusters and correspond- 
ing increase in average cluster size with humidity is 
consistent with the X-ray diffraction results described 
in the previous section and implies a cluster rearrange- 
ment similar to that proposed by Hsu & Gierke [16]. 
The percentage of the surface covered by clusters at 34 
and 9% humidities, is 17.6 and 18.9% respectively. The 
larger percentage area at lower humidities can easily be 
explained by tip broadening effects, whereby a greater 
number of small features can appear to have dispro- 
portionately larger area than fewer larger features. This 
is considerably larger than a value of 7.6% calculated 
assuming: that the membrane is homogeneous and dry, 
that the mass of SO3H is 8 1 . 1 g and the equivalent weight 
of the membrane to be 1070 (i.e. (81.1/1070) x 100). 
However, the membrane is not dry and will have incor- 
porated water preferentially into the hydrophilic clus- 
ters which will have swollen, the surface may not be 
representative of the bulk and the sulphonic acid group 
composition is 8% by mass, but not by volume. 



4. Conclusions 

The average number of water molecules per sulphonic 
acid group in Nafion 115 at 0%RH was found to be 
1.55. If it is assumed that all of the additional water is 
incorporated into the hydrophilic regions as the humid- 
ity increased then in the approximately linear region 
between (0-50) ±2% RH an increase of 14 ±2% in 
RH resulted in one additional water molecule per sul- 
phonic acid group. 

The most statistically likely structural model com- 
patible with the MaxEnt reconstruction of the small 
angle X-ray scattering (SAXS) data is an interparticle 
model in which point-like ionic clusters aggregate to 
form higher order agglomerates. 

Tapping mode phase imaging was successfully used 
to identify the hydrophobic and hydrophilic regions 
of Nafion perfluorosulphonate cation exchange mem- 
branes. Since there is often little correlation between 
the topography and phase images, it is a useful tool for 
identifying and mapping regions with different prop- 
erties, irrespective of their topographical nature. The 
images support the MaxEnt interpretion of a cluster 
model of ionic aggregation, with spacings between in- 
dividual clusters ranging from 3 to 5 nm, aggregating to 
form cluster agglomerates with sizes from 5 to 30 nm. 

The phase images also showed that the number of 
clusters decreased and the average cluster size increased 
with increasing humidity. This is consistent with the in- 
terpretation of the MaxEnt charge distributions, and is 
reminiscent of the redistribution of ionic material be- 
tween clusters proposed in the cluster-network model 
[16]. The change in number density of ionic clusters 
as a function of water content explains why the bulk 
volumetric swelling in water is observed to be signif- 
icantly less than the swelling inferred from scattering 
measurements. 



It is significant that the same conclusions were ob- 
tained using two different experimental techniques op- 
erating over complementary ranges of relative humid- 
ity. This indicates that the swelling and redisiribulion of 
ionic clusters in Nafion is a dynamic process occurring 
continuously from 0% to 100% relative humidity. 
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ABSTRACT 

Commercially available 4,4'-dichlorodiphenylsulfone (DCDPS) was successfully 
disulfonated with fuming sulfuric acid to yield 3,3'-disodiumsulfonyl-4,4'- 
dichlorodiphenylsulfone (SDCDPS). Subsequently, DCDPS and SDCDPS were 
systematically reacted with 4,4'-biphenol under nucleophilic step polymerization 
conditions to generate a series of high molecular weight, film-forming, ductile, ion 
conducting copolymers. These were converted to the acid form and investigated as 
proton exchange membranes for fuel cells. Hydrophilicity increased with the level of 
sulfonation. However, water sorption increased gradually until about 50 mole percent 
SDCDPS was incorporated, and thereafter showed a large increase to yield water soluble 
materials for the 100% SDCDPS system. Atomic force microscopy (AFM) confirmed 
that the morphology of the copolymers displayed continuity of the hydrophilic phase at 
60 mole percent SDCDPS. Conductivity measurements in the 40-50 mole percent 
SDCDPS range, where excellent mechanical strength was maintained, produced values of 
0.1 S/cm or higher which were comparable to the control, Nafion™. These compositions 
also show a high degree of compatibility with heteropolyacids such as phosphotungstic 
acid. These inorganic compounds provide a promising mechanism for obtaining 



conductivity at temperatures well above the boiling point of water and membrane 
compositions containing them are being actively pursued. 

The water soluble 100% SDCDPS system was further investigated by 
successfully functionalizing the endgroups to afford aromatic amines via appropriate 
endcapping with w-aminophenol. Oligomers and polymers from 5-30 kg/mole number 
average molecular weight were synthesized and well characterized by NMR 
spectroscopy, endgroup titrations and size exclusion chromatography. The diamino- 
telechelic sulfonated segment was reacted with several dianhydrides and diamines to 
produce muhiblock, hydrophobic polyimide-hydrophilic sulfonated polyarylene ether 
copolymers. Both ester-acid and amic acid synthesis routes were utilized in combination 
with spin-casting and bulk imidization. A series of tough, film-forming segmented 
copolymers was prepared and characterized. AFM measurements demonstrated the 
generation of quite well defined, nanophase-separated morphologies which were 
dependent upon composition as well as aging in a humid environment. Characterizations 
of the segmented copolymers for conductivity, and water and methanol sorption were 
performed and comparisons to state-of-the-art perfluorinated Nafion™ systems were 
made. It is concluded that the segmented or block systems have the potential to enhance 
certain desirable PEM characteristics in fuel cells, particularly those related to swelling, 
retention of mechanical strength at elevated temperatures, and critical adhesion issues in 
membrane electrode assemblies. 



films to remain stable during the chemical conversion from sulfonic sodium salt to 
sulfonic acid is crucial if the films are to be used as a PEM in a polymer electrolyte fuel 
cell. The triethyl ammonium salt form is needed to provide stability of the sulfonated 
monomers during polymerization, while the acid form is needed for proton conductivity 
when used as a PEM in a fuel cell. The polyimide containing 40 mole percent sulfonated 
diamine had an ion-exchange capacity very close to that of Nafion, while conductivity 
measurements yielded a value of 0.09 S/cm, which is near the value of Nafion at 0. 1 
S/cm. 

b. Chemical Modification of Polymers 

An alternate approach to introducing ionic functionalities into polymeric systems 
involves subsequent polymer modification reactions. This approach is especially suitable 
for the introduction of sulfonate groups via sulfonation of aromatic or unsaturated sites. 
The simplest example is the sulfonation of poly (styrene) by using either acetyl sulfate or 
sulfur trioxide/triethylsulfate, shown in figure 2.31.^°^ Acetyl sulfate is the reaction 
product of acetic anhydride and sulfuric acid. In the case of poly (styrene), a pendant 
sulfonate is formed, but main chain aromatic rings can also be functionalized with the 
sulfonic acid moiety. 

The sulfonation of ethylene-propylene-diene (EPDM) monomer systems has 
generated considerable commercial interest.'^"* By using Ziegler-Natta catalysis 
techniques, a random copolymerization of the three monomers can be achieved. One 
of the preferred dienes is 5-ethylidene-2-norbornene which contains an endo- and 
exocyclic group. Polymerization proceeds through the endocyclic double bond and the 
exocyclic unsaturated site is available for subsequent sulfonation, as shown in figure 
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Figure 2.3 1 : Sulfonation of Poly (styrene) 
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General Motors, Volkswagen, Nissan, Honda and Volvo among the largest participants in 
research and development. It has been projected that the number of vehicles (worldwide) 
will increase from the current level of 600 million to 1 billion by 2015.^^^ The American 
Methanol Institute recently estimated that by the year 2010, automakers will have 
introduced at least 2 million methanol fuel cell vehicles worldwide, with the number 
surpassing 35 million by 2020. Fuel cell vehicle introduction will focus initially on the 
three states in the U.S. requiring the sale of zero emission vehicles by 2003 (California, 
New York and Massachusetts) as well as Germany and Japan. These are areas of high 
population that are candidates for early acceptance of alternative fuel vehicles because 
they tend to have higher levels of pollution and offer maximum scale efficiencies for the 
first wave of methanol fuel infrastructure. 

The vast quantity of natural gas in the world ensures the availability of feedstock 
to produce the methanol needed for the future fleet of DMFC vehicles. In 1996, reserves 
stood at 49,912 trillion cubic feet (TCP) with annual consumption of 78 TCP. With a 
fleet of 35 million vehicles consuming 15.4 billion gallons of methanol annually (at 441 
gal/yr each), this would create a demand for 1.4 TCP of natural gas per annum, which is 
only 2 % of current annual consumption.*^^ It can be argued that the methanol supply 
and delivery infrastructure can be adjusted to meet future DMPC vehicle demands. 

7. Water Uptake and Protonic Conductivity in Sulfonated Polymer 
Membranes 

a. Water Uptake of Nafion Fluoropolymer 

Proton transport in ionic polymer membranes, especially when in the acid form, is 
largely influenced by the water content of the membrane. Ionic membranes behave like 
insulators in the dry state, but become conductive as a function of the water content when 
hydrated. Various studies have been conducted as to the minimum water content 
threshold for measurable conductivity. Several researchers have shown that a minimum 
of 6 or 7 water molecules are needed per sulfonic acid site to generate sufficient 

Lewis, R.; Doian, G. In Fuel Cell Power for Transportation; Seaba, J., Stobart, R., Eds.; Society of Automotive Engineers: 
Warrendale, PA. 1999; Vol. 1425, p 22. 
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conductivity for use as a PEM.^*^ Among the ionic polymers with high protonic 
conductivity, the activation energy of proton conduction of Nafion perfluoropolymer is 
low in comparison to other polymers* which may be due to the state of the water in the 
membrane. In a fuel cell using proton exchange membranes, the water content is 
controlled by humidification of the gas streams at each side of the membrane. Under 
these conditions, it has been shown* that the proton conductivity reaches a maximum at 
temperatures between 55-70*^C while water contents approach a minimum. This work 
established that at lower temperatures between 25-50°C, a variation in water content is 
less important than a temperature increase, which plays an important role in the kinetics 
of proton motion in the polymer membrane. 

Work by Zawodzinski et. al also involved investigation of water uptake of 
Nafion 117 membranes.*^^ Water uptake was studied by investigating membranes that 
were initially fully hydrated and were then exposed to a variety of drying conditions: 1) 
dried at room temperature under vacuum; 2) condition 1 followed by vacuum drying at 
105°C for an hour; 3) drying at room temperature over P2O5; and 4) condition 3 followed 
by vacuum drying at 1 05°C. Earlier work by Bunce et al has shown that the water 
content after drying under condition 1 resulted in X = 1 (where X = number of moles of 
water per mole of sulfonic acid group in polymer), whereas drying conditions 2, 3 and 4 
resulted in total dehydration. The extent of rehydration of the membranes was dependent 
on the drying condition. Figure 2.40 illustrates the water uptake as a function of time and 
compares a film that was dried at room temperature under vacuum versus a film dried 
under vacuum at lOS^'C. The membrane that was dried at room temperature contained the 
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Reike, P. C.; Vanderborgh, N. E. / Membrane Sci. 1987, 32, 313. 
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same amount of water after rehydration as it did prior to drying (X = 22), whereas the 
water uptake of the Nafion membrane dried at elevated temperature was only half of the 
initial value. Interestingly, the water uptake of the membrane dried at room temperature 
was independent of the water temperature, while the water content of the membrane dried 
at elevated temperature was a function of the water temperature. The results are 
summarized in table 2.2. 

The most likely explanation for the above results involves the formation and 
dissociation of the ionic clusters as the membrane is dried. The pretreatment of the 
Nafion 117 membranes were boiled in 0.5M sulfuric acid for an hour, followed by an 
hour in boiling deionized water to remove any excess acid. The pretreatment should 
provide enough thermal energy to form hydrated ionic clusters, but reorientation of the 
polymer during the drying step is required for the clusters to be broken up, which requires 
enough energy for Nafion ionic side chains motion. As long as water is present, the 
sulfonic acid moiety remains ionized, but as the membrane becomes dehydrated, the ionic 
clusters shrink in size and remain aggregated. As the final traces of water are removed, 
the sulfonic acid group is no longer ionized, the coulombic barrier to reorientation is no 
longer present, and the clusters can be dissociated at temperatures near the glass 
transition of the polymer. This explanation applies the same principles suggested by 
Eisenberg, where cluster formation is a balance between the ionic dipole interactions that 
favor formation and elastic forces that disrupt their formation.*^'* 

The dependence of membrane rehydration on drying conditions has important 
implications if they are used as an MEA in a fuel cell One method of MEA fabrication 
requires hot-pressing a pair of gas diffusion electrodes that containing a Nafion 
dispersion onto the surface of the Nafion 1 17 membrane at 120°C. This process causes 
complete dehydration of the membrane causing a potential decrease in X. This may have 
a negative effect on the maximum attainable conductivity of the membrane, since 
conductivity and water uptake have a roughly linear relationship. 
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b. Protonic Conductivity of Nafion Fluoropolymer 

Although many properties of Nafion fluoropolymer have been investigated*^^, 
proton conductivity details are reported less frequently. Some of the pivotal work has 
been performed by Zawodzinski et. al. '^*^and Kreuer et. al}^'^ 

The conductivity of Nafion 117*^^ has been reported by Zawodzinski et al. as a 
function of water content (figure 2.41) and temperature (figure 2.42). At 30°C, the 
conductivity decreases in a generally linear fashion with decreasing water content. The 
value of the conductivity, 0,06 S/cm, where ^ = 14 (where X = number of moles of water 
per mole of sulfonic acid group in polymer), agrees with data reported by Rieke eL al/^^ 
An Arrhenius plot illustrating the temperature dependence over the range of 25-90X is 
shown in figure 2.42. The observation that the plot is not linear illustrates the possibility 
that various mechanisms may be involved in protonic motion through the polymer 
membrane which can be calculated from the conductivity data. The diffusion coefficient 
of H"^ can be calculated from the conductivity by using the Nemst-Einstein Equation*"^^, 
below. This equation was derived by W. Nemst in 1888*^* and was originally used to 
calculate the diffusion coefficient of ions in dilute solutions. 
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Do = RT(r * r / AoXz"" + z ")(z'' * z ') 

Where: 



Do 


Diffusivity of molecule at infinite dilution (ionic pair), in cmVsec 


R 


Gas constant, 8.931 * 10"'*^ 


T 


Temperature, in Kelvin 




Cationic conductance at infinite dilution, in S/equivalent 




Anionic conductance at infinite dilution, in S/equivalent 


Ao 


Electrolyte conductance in infinite dilution, in S/equivalent 


+ 

z 


Valence of cation (absolute, no sign) 


2 ' 


Valence of anion (absolute, no sign) 
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A type of rewritable optical disk that employs the phase change recording method. Using 
this technique, the disk drive writes data with a laser that changes spots on the disk 
between amorphous and crystalline states. An optical head reads data by detecting the 
difference in reflected light from amorphous and crystalline spots. A medium-intensity 
pulse can then restore the original crystalline structure. 

Magneto-optical and dye-polymer technologies offer similar capabilities for developing re- 
writable optical disks . 
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Theories for the dielectric constant. 

G. van der Zwan 
December 4, 2003 

1 Maxwell equations and constitutive relations. 

The Maxwell equations in a polarizable, magnetizable medium, in which also charges are 
present, can be written in the following form: 

VB = 0 

V-^ = -(p-V-P) 

VxE + ^ = 0 
ot 

We neglect magnetization effects, thus setting M = 0. The following relation is always true: 

D = eoE + P (1.2) 

Therefore the equations can also be written as 

V-5 = 0 
V-^ = p 

WxE + ^ = 0 
at 

V.H-§ = r (1.3) 

where B = fioH in case the magnetization is zero. 

The charge conservation equation follows from taking the divergence of the last equation, 
using that V • (V x ^) = 0, and inserting the second equation. We thus get 

| = -v.T . (1.4, 

In the remainder of this section we use Fourier transform language, meaning that all partial 
derivatives with respect to time are replaced by —iuj. 
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We now assume the following two constitutive, linear relations: 

j{f,uj)=^K{u;)Eif,u) (1.5) 

and 

D{f,u;) = eoe{uj)E{f,uj) (1.6) 
Eliminating density /?, current j, and D in virtue of E, we can write these equations as 

V X E{f,u) - iujfioH{f,uj) = 0 (1.7) 

and 

V X H{f,u;)'\-iuJ (^eoe{u) - E{f,uj) = 0 (1.8) 

whereas for the charge density we find that 

p(f,u) = '^^-E{f,u) (1.9) 

Taking once more the curl of eq. (1.8), using that 

V X (V X ^) = V(V 'H)-AH= -AH (1.10) 
in view of the first of eqs. (1.3) and the absence of magnetization, we find: 

-AH{f,u)+iuj (^eoe{uj) - V x E{f,uj) = 0 (1.11) 

Subsequent substitution of (1.7) the yields: 

AHif.uj) + ^ - ^\ H{f,u) - 0 (1.12) 



This is a wave equation for the vector field with dispersion relation: 



We can derive a similar equation for the electric field: take the curl of eq. (1.7) and use the 
equivalent of relation (1.10), we get: 

AE{f,u;) + ^ (e{u) - ^) E(r» = -^Vp(f,c.) (1.14) 

where we also used eqs. (1.8) and (1.9). 

This is again a wave equation, now for the electric field with the same dispersion relation 
as above, but with a source term. 

The assumption is now that the charges move slow compared to changes in the electro- 
magnetic field. In that case we can neglect the source term. Alternatively we can state 
that the transverse waves are unaffected by the charge distribution, the longitudinal wave 
probably decay faster as a function opf the distance to the sample. A bigger problem is 
the frequency dependence of k. Assuming that it can be replace by its zero-frequency 
value means that the charges can follow the electric field. For microwave frequencies that 
is 10^^ — 10^^5"^ this seems strange. 
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2 The dielectric constant of a gas of non-polarizable 
molecules. 

The simplest theory for the dielectric constant is that for a gas consisting of non-polarizable 
dipoles. 

The dielectric constant can be related to the ratio of the polarization, or macroscopic dipole 
moment per unit volume of a system, and the electric field E in the following way: 

P{r) - eoxm = eo(6 - l)Eif) (2.1) 

Note that P is dipole moment per volume, which makes x ^ dimensionless quantity. 
In reahty x gives the relation between electric field and polarization, where the electric field 
E is not necessarily equal to the external field Ee, but is modified by the presence of a 
polarizable medium. 

We therefore have to calculate the polarization for a given external field to obtain e, as well 
as the internal electric field. 

For a low density gas of polar, non-polarizable molecules the interaction between the 
molecules can be neglected, and we only need the response of one dipole to an electric 
field. 

For the moment we assume that the internal field is equal to the external field. Since we 
are in an equilibrium situation, the probabihty of finding a dipole ft with given orientation 
in a field E^ is determined by the Boltzmann factor: 

Q^flEe COS l9 

where t? is the angle between dipole and field. The average dipole moment is then determined 
by 

(cost?) = 1^ cos7?cosz9P(cosi9) = (2.3) 

where we used the small field approximation, that is, expanded the exponential to lowest 
relevant order. 

This means that we can write the average dipole moment as: 

(u) - -^^e (2.4) 
Consequently the polarization of the gas is given by: 

where p is here the density of the gas. Therefore the dielectric constant is given by: 
again assuming that internal and external field are equal. 



One wonders whether this can be taken a step further. After all the polarization causes an 
electric field, so that the field actually felt by the dipoles is_diflFerent. If there is a polarization 
P, every volume element dV = df has a dipole moment Pdf. A dipole at position r causes 
an electric field at position f ' given by: 

^ '^l - 3^ V P(r-)dr (2.7) 



where 

R = f~f' (2.8) 

Since the system is infinite, we can calculate the field at the origin, and since the polarization 
is constant, the integral can easily be calculated: the field due to polarization is equal to 



Ep — 
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P' jdn(l~ SnO) jdr^ (2.9) 



This does not work very well, unless we exclude a sphere or something around the origin to 
avoid infinities. The reason is that in this method we allow polarization at the position we 
measure the field. Anther possibility would be to take individual dipoles, connected to the 
molecules, and perform an average over the positions of the molecules. This also leads to 
infinities, since there is always a finite probability of finding the molecule at the origin. In 
addition we have a problem for large values of r, which is a consequence of the long range 
character of the electrical forces. No easy way to get around that either unless we want to 
take the shape of the vessel in which this is all taking place into account. SO: problems all 
around with this approach. See in this context also the remark by Fixman [1], p. 2076. 
Let's try the probability method though, and see what happens. Before that the following: 
FYohlich [2] assumes in the derivation of the gas phase dielectric constant that the dipoles do 
not interact, on the bases of the assumption that if they are far enough apart, the interaction 
energy is small compared to ksT. This means that 

« ksT (2.10) 



47reor' 

where r is the distance between the dipoles. At room temperature ~ 298 K, we have fc^T ^ 
4.11 X 10-23 J. For a 1L> - 3.336 x 10"^^ Cm dipole we get fi'^/Aireor^ ^ lO'^Vr^ Jml 
This means that for the above condition to hold we must have: 

r » 1.4 X 10"^ m - L4nm (2.11) 

In fact we want the electric field due to the dipoles at the position of the central dipole also 
to be small compared to the external electric field. 
The microscopic dipole density is given by: 

p(r)=f:M^"'-^"'i) (212) 

1=1 

where pi is the dipole moment of particle i and u its position. 
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What we are interested in is the electric field at the position of a dipole due to all the other 
dipoles. The electric field at position f due to a dipole at position f\ is given by: 

m- , -,3 fi-3 'V-":,r'' Vfi (2.:3) 

47r6o|r -ri|3 \r - Ti\^ J 

Thus the total field at the position of dipole i is given by: 

The next step is to perform averages: both over all positions of the dipoles, and over 
their orientation. Since we neglected interactions between the dipoles, they are uniformely 
distributed around the central dipole. 

Using result (2.4), but now with the orienting field equal to the sum of field due to the 
polarization and external field, we get 

where we also used that the magnitude of all dipoles is equal. The remaining average is 
over the positions of the dipoles. This average is therefore defined as: 

{■■■) = jdr^-'P{{r^-'}){---) (2.16) 

where P{{f^~^}) is the probability of finding a configuration {r^~^}. Note that strictly 
speaking we need the conditional probabihty of finding a configuration of — 1 particles, 
given a dipole at position fi, but since we neglected interactions between the molecules 
these are all independent, and in fact the probability of finding a particle anywhere is just 
where V is the volume of the system. Thus we can write the average as: 

To get this result we used that all particles are equivalent, and therefore each of the terms, 
of which there are N — 1, in the sum is the same, and furthermore that all positions are 
equivalent so that we can place the central particle to the origin. 

In fact looking again at the earlier result, and this result, we notice that the angular inte- 
gration equals zero. This can be understood in the following way: the probaihty of finding 
a dipole with given orientation at position fi is equal to finding a dipole with the same 
orientation at position — fi, since the orientation is assumed to be uncorrelated with the 
orientation of the central dipole. The field due to these dipoles is equal, but opposite, and 
therefore the total contribution to the field is zero. The whole exercise therefore showed 
taht we can indeed take the external field here as long as we neglect correlations between 
the dipoles. 

But now I have a problem. This calculation of the inside field shows that there is no 
contribution due to the polarization when the field is constant. But it is always stated that 
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the electric field in a dielectric is reduced by the dielectric constant, so we also have the 
following for the field between two charged plates with charge density a. The external, or 
vacuum field, is given by 

E, = - (2.18) 

whereas the field in the dielectric is 

E = — (2.19) 

This means that the polarization must give rise to the difference between the two fields: 

Ep^E-E^ = — ^ (2.20) 

So according to this the electric field does change due to the presence of dipoles. If it is a 
correlation effect, than why are we allowed to neglect it in the derivation? 
Another remaining problem, also present in the last derivation, is that, if the angular inte- 
gration does not make the whole polarization effect vanish, we are still stuck with an infinite 
contribution due to the long range nature of the electrical forces. 



3 The dielectric constant of a gas of polarizable molecules. 

The first refinement can be made by making the molecules polarizable. This means that 
the dipole itself responds to the field by becoming larger. This enlargement is measured by 
the polarizability a. In that case the average dipole moment in an external field is given by 

Thus we get for the static dielectric constant: 



Introducing the high frequency dielectric constant eoo as 



- 1 = ^ (3.3) 

we can write the previous expression in the form: 

4 A dilute solution of dipoles in a polarizable medium. 

The next step is to consider a dilute solution of dipoles in a polarizable medium. The 
dipoles are so far apart that they still do not have direct interaction, but since the medium 
is polarizable, we have to account for the fact that the internal field is different from the 
external field, and in addition that the reaction field of the dipoles influences the dipole 
moment if they are themselves polarizable. 
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5 Dielectric constant according to Debye. 

Debye [3] starts out to find an explanation for the temperature dependence of the dielectric 
constant. Up to that time it was apparently customary to think of the dielectric behavior as 
resulting from the displacement of charges within an isolator, and then the dielectric constant 
does not depend on temperature. Rather than abolishing statistical mechanics Debye prefers 
to assume that in an isolator we do not just have bound charges but also permament dipole 
moments, and with that hypothesis he can explain the temperature dependence, on the 
basis of classical statistical mechanics. 

On the basis of that hypothesis he basically derives eq. (3.2) with one important difference: 
for the field actually felt by the molecule, the directing field he does not use the external 
field, but a field derived by Lorentz (the so-called Lorentz field) which is given by: 

E = E, + —P^^E, (5.1) 

which then, with the same calculation for the avearge orientation of the dipoles leads to the 
expression: 

e + 2 Qeo/CBT 3eo ^ * ^ 

or, with the introduction of the high frequency dielctric constant, for which a similar argu- 
ment can be used: 

e + 2 eoo + 2 ^e^ksT ^ ' ^ 

Debye then continues to show that for a number of substances the graph of the function 
(e - l)r/(6 + 2) vs T is indeed a straight fine. 

The internal field according to Lorentz. 

Lorentz [4] gives a derivation of the internal field which we will quote entirely adapting the 
notation at various places, as to conform with the notation introduced here. 

We shall observe in the first place that the field in the immediate neighborhood 
of a polarized particle may be determined by the laws of electrostatics, even 
when the electric moment is not constant. Take, for instance, the case treated 
in § 43 (this paragraph treats dipole radiation). It was stated in note 23 that 
at great distances the terms resulting from the differentiation of the goniometric 
function are very much greater than those which arise from the differentiation 
of K These latter, on the contrary, predominate when we confine ourselves to 
distances that are very small compared to the wavelength; then we may write 
(for the potentials of the time dependent dipole p): 

^ -p- and a = 0 (1) 



47reo T 

from which it appears that the field is identical with the electrostatic field that 
would exist, if the moment p were kept constant. 
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It is further to be noted that the difference between the mean electric force and 
the electric force existing in a small cavity depends only on the actions going on 
at very small distances, so that we may deal with this difference as if we had to 
do with an electrostic system. 

Let us therefore consider a system of molecules with invariable electric moments 
and go into some details concerning the electric force existing in it. 

The field produced by the electrons being determined by 

A0=--^, and e=-V0 (2) 
we have for the mean values 

A0--^, and E = ?--V^ (3) 
Co 

or, in words: the mean electric force is equal to the force that would be produced 
by a charge distributed with the mean or, let us say, the "effective" density p. 

In the definition of a mean value 0 given in §113 it was expressly stated that the 
space S was of spherical form. It is easily seen, however, that we may as well 
give it any shape we like, provided it be infinitely small in a physical sense. The 
equation 

pS- J pdS (4) 

may therefore be interpreted by saying that for any space of the said kind the 
effective charge (meaning by these words the product of (p and S) is equal to the 
total real charge. 

We shall now examine the distribution of the effective charge. Suppose, for the 
sake of simplicity, that a molecule contains two electrons situated at the points 
A and B with charges — e and +e, (?) and denote by r the vector AB. There 
will be as many of these vectors, of different directions and lengths, as there 
are molecules. Now, if the length of these vectors is much greater than the size 
of the electrons, we may neglect the intersections of the bounding surfaces of 
the space S with the electrons themselves, but there will be a great number of 
intersections with the lines AB. These may not be left out of account, because 
for any complete molecule / pdS = 0, whereas each of the said intersections 
contributes to the effective charge within a (the bounding surface of S) by an 
amount -e or +e as r*n (where n is the normal to a drawn outwards) is positive 
or negative. Hence the total charge within a may be represented by a surface 
integral. In order to find the part of it corresponding to an element da (infinitely 
small in a physical sense) we begin by fixing our attention on those among the 
lines AB which have some definite direction and some definite length. If the 
starting points A are iregularly distributed and if, for the group considerd, there 
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number per unit volume is ly, the number of interseections with da will be ur-nda 
when r • n is positive, and —yr- ftda when it is negative. Therefore, the part 
contributed to the charge within a is —i/ef-nda in both cases, and the total part 
associated with da is —T^uf- ftda, the sum being extended to all the groups 
of lines AB. But ef is the electric moment of a particle, uef the moment per 
unit of volume of the chosen group, and J2 the total moment per unit of 
volume. Denoting this by P, we have for the above expression — ^ vef • n the 
value —P'fi, and for the effective charge enclosed by the surface a 



J 



P-fida (5) 



As the difference between E and the electric force in a cavity depends exclusively 
on the state of the system in the immediate vicinity of the point considered, we 
may now conceive the polarization P to be uniform. In this case the above 
integral is zero for any closed surface lying interely within the body, so that the 
effective charge may be said to have its seat on the bounding surface S. Its 
surface density is found by calculating the above integral for the surface of a 
flat cylinder, the two plane sides of which are on both sides of an element dL 
at a distance from each other that is infinitely small in comparison with the 
dimensions of dS. Calling N the normal to the surface E, we have at the outer 
plane P • n = 0 ( if we suppose the body to be surrounded by ether), and at 
the inner one P • ft = -P • N. The amount of effective charge contained in 
the cylinder is therefore given by P • TVdS, and the charge may be said to be 
distributed over the surface with a density P • N. 

Now consider a point A of the body. By what has been said, the electric force E 
at this point is due to the charge P - N on the bounding surface S. If, however, 
a spherical cavity is made around A as center, there will be at this point an 
additional electrioc force E\ caused by a similar charge on the walls of the 
cavity, and obviously having the direction of P. The magnitude of this force is 
found as follows. Let a be the radius of the sphere, da an element of its surface, 
1? the angle between the radius drawn towards this element and the polarization 
P. The surface density on da being -Pcost?, we have for the force produced at 
A 



giving 



^ / Pcos^dda (6) 



E' = \P (7) 



Our foregoing remarks show that the expression 



E + \P (8) 
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may always be used for the electric force at the center of a spherical cavity, even 
though the polarization of the body change from point to point and from one 
instant to the next. 

The reasoning for getting eq. (6) is rather unclear to me. 

The previous quote shows that Lorentz attributes the field in the cavity to the surface 
charge on a spherical cavity. Let's try a diflFerent derivation. Suppose the polarization in 
the medium is given by P{f). The potential at any point is then given by the integral 
equation: 

This integral is over all space, and in view of the discontinuity in P at the surface of the 
cavity, the divergence has a delta function, and we have to be careful evaluating this integral. 
This can be converted to a surface integral in the following way: first we write: 

We can now use Gauss' theorem to rewrite the first term as a surface integral: 

^^r> = ^ / - V . / ,f'p£l (5.6) 

Note that the volume is outside the cavity, so the normal points "inwards" . 

In the second integral there is now no problem taking P as a constant, since that does no 

longer lead to singularities. 

Putting P a constant everywhere then yields: 

0(f) = — P • / ds'—^ - — P • V / (5.7) 
47reo Js \r — r'\ Airco Jv |r — r'| 

We want to caluculate the electric field, so we take the gradient, and we want it at the center 
of the cavity, so we put r = 0, but after taking the gradients. This leads to the expression: 

E=- (VcPif)) ^ = -r^P • / drn'n' + -^P • / df 'V'V'^ (5.8) 
V Vr=o 47reoa2 Js 47reo Jv r' 

where the second integral gives 0, for the same reasons as were explained in section 2. The 
first integral gives ^o?, so that 

E, = ±.P (5.9) 

is the electric field at the origin of the cavity due to the polarization. This is exactly Lorentz' 
result. The implication is that the effect of the boundary is partially taken into account. In 
fact one should solve the boundary value problem first, and then calculate the internal field 
with the above integral. 

In this context: see also Lorentz[4], §117, and Prohlich[2], p. 22. 
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6 Dielectric constant according to Onsager. 



Next we come to Onsager's treatment of the dielectric constant. Onsager [5] goes into 
considerable more detail, which will be outlined in this section. 

His major critique is that Debye uses the total internal field to orient the molecule, whereas 
according to his way of thinking about it, only part of this field should be used, namely that 
part which does not depend on the original dipole itself (and which of course contributes to 
the polarization). So he needs to calculate two fields: the field due to the dielectric outside 
the cavity, and the reaction field due to the dipole inside the cavity. His starting point is 
that the Clausius-Mosotti (or Debye) expression eq. (6.2) leads to a phase transition for the 
polarization, namely for a certain critical temperature Tc, e goes to oo, which is only very 
rarely obeserved, but it also means that fluids with high dielectric constants are necessarily 
close to this phase transition, and should therefore exhibit large (critical) fluctuations. This 
is also not observed. 

The first step is to derive the reaction field of a dipole at the center of a spherical cavity, 
again with radius a. 

The solutions for the inside- and outside electric fields due to such a dipole ft can be written 
as 

E,^-^p-{l-3fr)-R (6.1) 
where the reaction field R is given by: 

--^^'-^^ <-) 

and 

Eout - T^f^ ■ (1 - 3rr) (6.3) 

where the external moment /T is given by 

? = ^.- (6.4) 

The second boundary value problem that has to be solved is the one alluded to at the end 
of the previous section, namely that of an empty cavity, embedded in a dielectric in which 
a constant electric field E is present. This boundary value problem leads to the following 
expressions for internal and external fields: 

E^^G = -^E (6.5) 
2e + 1 

and 

4ut = • (1 - 3rf) + E (6.6) 
2e + 1 H 

The total field acting on the dipole in a spherical cavity is therefore 

F = G + R^-^E+ ^ n (6.7) 

2e + l 47r€oa^(2e + 1) 
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and the dipole moment of the molecule is given by 

il = po + aF (6.8) 

where jh is the permament dipole moment of the molecule in the absence of external fields. 
Furthermore we can introduce the refractive index n through the relation 

a 3 n2-l n2-l 
« (6.9) 



47reoo^ ^na^p 'n? + 2 v? -\-2 

With this relation it is easy to derive from eqs. (7.7) and (7.8) the dipole moment of the 
molecule: 

^ (2e + l)(n2 + 2)„ e{n^-\), 

The next step in Onsager's calculation is establishing the potential of mean torque. He gives 
two methods, of which we will only give the first. The basic problem is that the reaction 
field should not contribute to the torque on the dipole since it is in the direction of the 
dipole moment. The situation is made complicated by the polarizabihty. But it is fairly 
simple to calculate the torque on the dipole fi, which is of course given by: 

M = Fxp (6.11) 

and subsequently find the potential which will give rise to this torque. (NB. compare to free 

energy calculations). 

A direct calculation shows that 

FxiI=-^^Exi2 = Exfr (6.12) 
2e + l ^ ^ 

We note, however, that the part of fi proportional to E does not contribute to the torque, 
and by inserting (7.10) we obtain: 

SO that 

^ e(n2 + 2)^ 

^° = -277^^^ ^^-''^ 

This means that the potential of mean torque is given by 

w = ~fro^E (6.15) 

which is the potential that has to be substituted in the Boltzmann factor. We note here 
that, if the polarizabihty is neglected, this is just the internal field, but now solved with the 
correct boundary condition, whereas Lorentz did not bother to solve the boundary value 
problem, just took the polarization constant everywhere. 
This allows us to calculate (il-E). Using (7.10) we get: 

(fl-E)=^-^ {it, . E) + '.^^4neoa^E (6.16) 
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and with (7.15) it is easy to see that (see also the calculation leading to (2.3)): 




(6.17) 



Using this we find for the polarization in the medium the following expression: 



/ e(n^ + 2)^(26 + 1) nl 6(n^ + 2) \ ^ 
V (2e + n2)2 QfeeT 2t + v? j 



(6.18) 



We therefore find for the dielectric constant the following consistency relation: 




e(n^ + 2)^(2e + l) /zg e(n^ + 2) 
(2e + n2)2 9A;bT 2e + v? 



a 



(6.19) 



Again a strightforward calculation then yields the following relation: 

{2e + n^){e-v?) _ p^l 



(6.20) 



e{v? + 2)2 Qeofc^T 



It is to be noted here that for the polarizabihty part, the Clausius-Mosotti relation is again 
used, which seems strange, since the internal field for the derivation of Clausius-Mosotti 
is the Lorentz cavity field, and not the internal field found by solving the boundary value 
problem. 

This ends our exposure of Onsager's theory. He has a few more things to say, especially 
about mixtures, but these are not of immediate interest to us here. At the end of this 
section we also note that so far all theories are macroscopic, and basically try to derive 
expressions for the dielectic constant using consistency arguments. This will change in the 
next sections, where we start out from fundamental statistical mechanics considerations. 

7 Dielectric constant according to Kirkwood. 

In this section we outline the theory of dielectric polarization of Kirkwood [6]. Kirkwood's 
generalization is an extension of Onsager to take better into account the local correlations. 
A quote from Kirkwood reads: 

The most serious defect of the model lies in the assumption of a uniform lo- 
cal dielectric constant, identical with the macroscopic dielectric constant of the 
medium. Moreover, the hindered rotation of the neighbors of a molecule, rel- 
ative to itself, which is implicit in the deviation of the local from the average 
polarization, is attributed entirely to the dipole field of'the central molecule. 

We read this to mean that the reaction field, as calculated by Onsager is inadequate, since 
in the immediate neighborhood of the central molecule (cavity), due to for instance large 
electric fields the dielectric behavior is not equal to the bulk dielectric behavior, and secondly 
that short range molecular forces can also be of importance in determining the average 
orientation of a molecule close to the central molecule. 
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Starting point of Kirkwood's derivation is the classical statistical mechanical equihbrium 
partition function of a collection of dipoles in an external field. The total electric moment 
M is the vector sum of the molecular dipole moments jli: 

M^Zili (7.1) 

and the average dipole moment component of molecule i in the direction of the field E, 
denoted hy (^pi • E'j is defined as: 

where V/v is the total A^'-particle interaction potential (including short range and dipolar 
interactions), and dT a phase space element. 

Expansion in powers of the electric field then yields to lowest order: 

{il,.E)^ = 0{{il,-E){M-E))Eo (7.3) 

where the average on the right hand side is now to be calculated in the absence of of an 
external field. 

Next Kirkwood assumes that the internal field can be related to the external field by the 
relation 

E = (7.4) 

which means that he takes a spherical dielectric medium and shoves it into the external 
field; secondly, using isotropy of the system in the absence of external field, we can write 

The remaining average is performed in two steps: first over all configurations for given fli^ 
and subsequently over all orientations of fii. Since the averages are performed in the absence 
of fields, this last step is irrelevant. Thus we find: 

{fi,-M)^{jii-M^^il-^ (7.6) 

In this expression Mi is the average of M for given fli^ and in the final step we used that 
this average must be the same for all molecules, so the index can be dropped. 
This all leads to the Clausius-Mosotti type expression 

^ £^ (7.7) 
€ + 2 9eokBT ^ / 

which indeed reduces to Clausius-Mosotti if for M we substitute /i, in other words if we take 
all dipoles to be uncorrelated, except for the central dipole itself of course. Also Onsager 
should result if the correlation is given by just the reaction field. 
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Using the results of appendix A, we derive for the correlation: 

— 9e 2 



ii-M = -— (7.8) 

(e + 2)(2e + l) ^ ' 



Introducing this into eq. (8.7) yields: 

(2e + l)(e-l) p^i^ 



eoksT 



(7.9) 



which is indeed eq. (7.20) with = 1. 

Kirkwood goes one step beyond this result, by taking local correlation effects into account. 
The above derivation of Onsager's result shows that Onsager takes at least part of this 
correlation into account, namely the reaction field, but according to Kirkwood this is not 
sufficient, (see the above quote) . 

Kirkwood's improvement consists of the fact that rather than directly using eq. (A. 8) he 
states that this result is in fact only valid if no short range contributions to the interaction 
potential, and no saturation effects are important. 
He derives the following expression: 

(2e + l)(6-l) ^ pfl'p 

e eoksT ^ ' ^ 

where fl now contains the effect of short range correlations. An explicit expression is given 
for the case where there is a rather welh defined rigid solvation shell around the central 
dipole: 

/2-^ = /i2(l + ^{cos7)) (7.11) 
where the solvation number is z and (cos 7) is the average orientation between two dipoles. 



8 Dielectric constant according to Cole. 

One thing missing in Kirkwood's derivation is the effect of polarizabihty, although at various 
places (for instance the calculation of the central dipole moment) it is taken into account on 
an ad hoc basis. Cole [7] tries to improve on that on the basis of a simple harmonic oscillator 
model. Apparently a number of attempts were made to include polarizabihty directly, 
leading to different expresions, but they missed the consistency of Onsager's approach. (See 
other entries in te literature list). The main problem appears to be that calculation of the 
polarization for given central dipole on the basis of statistical mechanics leads to a conceptual 
problem, in that the dipole moment of a polarizable dipole cannot be fixed without fixing the 
electric field, meaning the dipole moments of all the other dipoles, and this is not necessarily 
an equihbrium configuration. Cole gets around this problem by explicitly taking the internal 
coordinate into account, albeit in the harmonic oscillator approximation. 
So, Cole's starting point is the interaction Hamiltonian: 

^ = Us-^lT.f'i^i+l E(/^oi + eifi) . Tij • {poj + ejTj) - ^(/Zoi + eifi) • Eo (8.1) 
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where Us is a short range interaction energy between the molecules which is not further 
specified, fi is an internal coordinate, which oscillates harmonically with spring constant ki] 
/?0z is the permanent dipole moment of molecule i, and Cifi the induced moment, so that the 
third term represents the dipole-dipole interaction, Tij being the dipole tensor, and the last 
term is the dipole interaction energy with an external field. Finally the tensor T is given by 

where Rij denotes the connecting vector between the particles i and j. 
The force on the coordinate fi is for a given field Eq given by 

kifi + eiY^T,j (8,3) 

where we only took the dipolar part due to the other polarizabilities, not due to the per- 
manent dipoles. Supposedly this is done in van Vleck, but so far I have not seen that. 
The equilibrium value for the r^, f^ij in the presence of an external field is given by 



9 Final remarks. 

1. This is something I don't understand: take the refractive index in the Onsager deriva- 
tion equal to 1. In that case the potential of mean force: eq. (7.15) is given by: 

w = --r^^Po 'E^ -flo' Ein (9.1) 

and the reason is that the internal field is now calculated from a boundary value 
problem, cf. eq. (7.5). 

This implies that the desired average is given by: 

/p.E) = — f^E (9.2) 

so that the polarization becomes: 

P^Pld^E (9.3) 



(2e + l)kBT 
which leads to 

(2e + l)(e-l) piil 



(9.4) 



which is again eq. (7.20) with — \, 

The point we can make here is the following: it is only the internal field that is changed 
from the Lorentz field, |(e + 2)E to the internal cavity field ^^E. No reaction field 
part enters here. 
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On the other hand, however, we have shown that Kirkwood's method, combined with 
the dipolar reaction field as the correlation between the central dipole and the other 
dipoles leads to exactly the same result. Here we do not have to calculate the internal 
field, only correlation functions in the absence of a field. Somehow this correlation 
appears to be missing from Onsager's work. 

It could be a coincidence based on a relation Onsager points out between on the one 
hand the internal field G and eq. (7.5), which yields 

- = — (9.5) 
E 2e + l ^ ^ 

and on the other hand the relation between external dipole moment and vacuum dipole 
moment 

1^ = ^ (9.6) 
fi 2e-\-l ^ ^ 

which are exactly the same. 

(Side remark: maybe worth looking into in the case of quadrupolar dielectrics?). 

In general Clausius-Mosotti is considered a good expression for gases, or for polarizable 
media with no permanent dipole moments on the molecules. Nevertheless it is not true 
that if we let the density in Onsager's expression, which, with Kirkwoods modification, 
is considered good for polar media, become small, we get Clausius-Mosotti back. This 
fact is probably related to something else: if we solve the boundary value problems, 
and finally let the cavity radius go to zero, we get a different result than for the case 
where we begin with no cavity at all. If we have no cavity at all, the internal field can 
be written as: 

4n = -E (9.7) 

e 

which leads to: 

an expression I have never encountered. 

In addition, for polarizable media it could be argued that the cavity is filled with the 
same polarizable molecule that the fluid consists of, so that (10.7) would indeed be a 
better expression for the internal field. For polarizable media we would then get 

eoo(eoo - 1) = ^ (9.9) 

^0 

In this context, van Vleck [8] claims the following 

Our method of partition functions is to be contrasted with the usual, es- 
sentially static Lorentz method of representing dipole-dipole coupling by a 
local field (• • •) for a long test body. The Lorentz procedure is shown to be 
only a first approximation, which is really warranted if the density is so low, 
or the temperature so high that one may neglect all terms but the first in 
the development of the partition function in 1/T. 
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The question arises if the boundary conditions need modification in the case of low 
density, and if so, how? Look into the development of hydrodynamic boundary con- 
ditions? 

See in this context also appendix C. 
Consistency. 

Another thing that is slightly mysterious is the following: To calculate the relation 
between polarizability, the simplest model is just to take the induced dipole moment 
proportional to the vacuum field, and the susceptibility the relation between the po- 
larization and the vacuum field: for a gas of polariable non-dipolar molecules we then 
get (cf. eq. (3.3)) 

eoo - 1 = ^ (9.10) 

^0 

This seems consistent. The next level of refinement is to take the Lorentz field as 
the field acting on the molecule, thus inducing polarization, which leads to Clausius- 
Mosotti: (cf. eq. (5.3)) 

This seems somehow inconsistent, in that the susceptibility is now taken to be the 
proportionality factor with the vacuum field, whereas it should be the prportionality 
factor to the internal field. If we do that, however, of course eq. (9.10) would be 
recovered. 

Now the next level of refinement would be to use the Onsager cavity field for the 
internal field. In that case we would find, by the same (inconstent?) procedure: 

(26oo + l)(eoo-l) ^ poa 

which is not what Onsager uses for the high polarizabiUty, in stead he uses Clausius- 
Mosotti, which seems even more inconsistent. 

The results of appendix C, however, show the correctness of Clausius -Mosotti for low 
densities, or whatever the exact expansion parameter is in that procedure. 

We can also look at the expansion of eo© as a function of x = pooc/eo for the various 
cases: Eq. (9.10) obviously leads to: 

e^ = l+x (9.13) 
Expansion of the Clausius-Mosotti function eq. (9.11) gives: 

eZ^ = l+x + lx' + ^x'-.. (9.14) 

and finally Onsager (9.12) also yields 

eZ"' = l+x + ^x'-^x' + --- (9.15) 
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which shows that Onsager and Clausius-Mosotti agree up to second order in the 
expansion parameter, 

4, Correlations. 

In the Lorentz model, there are no correlations between the existence of a hole and 
the polarization field outside. Still there is an internal field diflPerent from the external 
field. In the Onsager model, the presence of the hole causes the polarization to be 
modified. The hole "creates " a dipolar polarization field around it, and this dipolar 
field again leads to the modification of the inside field. We can therefore view the 
difference between the Lorentz and the Onsager field inside as a correlation eflFect: 



It seems to me that this must be sort of a maximum. For low densities, the system 
does not even realize there is a hole, the distances between the molecules are such that 
holes the size of the particle radius are constantly present. The other extreme is that 
the distances between the molecules are so small that all space is filled with dielctric. 
So why doesn't the density enter into this part of the picture? 

Now here is something I had not realized before. Lorentz also neglects the correlations 
between the existence of a central dipole and the outside polarization. Onsager puts 
a polarizable dipole in the cavity, this dipole polarizes the surrounding and therefore 
causes a reaction field which again acts on the dipole, thus enlarging it. Calculation of 
this effect, gives an additional contribution to the internal electric field. This can be 
seen as follows: The dipole moment is proportional to the total internal electric field: 




8e-2 



E 



(9.16) 



3(2e + 1) 




(9.17) 



The reaction field R is given by eq. (6.2): 

^ _ 1 2(e - 1) ^ 
~ 47reoa3 2e + 1 ^ 

Combining these last two equations the gives: 



(9.18) 




(9.19) 



or, 




(9.20) 



and hence 




(9.21) 



If we now assume that the density is related to the radius by 




(9.22) 
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and using that P = pa = eo(e - 1)E we then find the following relation 



3epa 



= eo{e - 1) 



(9.23) 



2e + l-|ir(^-l) 



from which we again find Clausius-Mosotti. So the remark made under the previuos 
note that there is an additional inconsitency in Onsager is in fact incorrect. 

Yet another way of looking at this is that the combined effect of the hole-polarization 
correlation, and the reaction field business again yields the Lorentz field [2]. Using 
Clausius-Mosotti to eliminate pa in (9.21) we find after some simple calculations: 



which is indeed the Lorentz field. 

The areas of application are different: Lorentz is supposed to hold for low density 
gases, whereas for Onsager (9.22) has to be used, which is in fact a condition for all 
space to be completed filled with spheres. Another coincidence with which this whole 
field abounds. 



(9.24) 
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A Dipole in a cavity at the centre of a dielectric sphere. 

Consider the situation where we have a spherical cavity at the centre of a dielectric sphere. 
The radius of the cavity is assumed to be a, and of the sphere it is given by R. We have 
potentials in three different regions, given respectively by 

ipi = — — T + Au'f, for r <a (A. 25) 



in the core region, 

fl • r 

in the intermediate region, and finally 



iP2 = B-P — -+Cp'f, for a<r<R (A,26) 



^3 = D-f^, for r>R (A.27) 

The boundary conditions determine the values of the four coefficients. The relations between 
the coefficients are: 

Aneoa^A- B - Aneoa^C = -1 
47reoaM + 2e5 - Aireoea^C = 2 
B + AireoR^C-D = 0 
-2e5 + 47reoei?^C + 2D = 0 (A,28) 

The solutions for the coefficients are: 

Aneoa^A 
B 



2(l-e) + 26|,f^ 
2^ + 1-2^^ 



6 + 2 

D - —5 (A.29) 

To calculate the necessary integral we need the polarization in the dielectric medium, which 
is given by: 

P - -6o(e - 1)VV'2 - - ^^y^ P • (1 - 3fr) - eo(e - 1)01 (A.30) 

In the integration over the medium the first term again cancels in view of the symmetry, 
and we find: 

drP = -6o(e - l)/xC J^^^ df = - yeole - mCR" f 1 - j (A.31) 

21 



In the limit R —* oo this becomes: 

hm f drP = {A.32) 
R^ooJ (e + 2)(2e + 1)'^ ^ ^ 

So with this approximation M becomes: 

At the end of this section we remark that it was necessary to first consider a finite region 
and finally take the limit to an infinite medium to avoid the infinities alluded to in previous 
sections. 



B Kirkwood's derivation. 

In this appendix we paraphrase Kirkwood's derivation of the reaction field, to show more 
clearly where the differences with Onsager are. For the larger part the derivation in the 
previous appendix can be used. The major change is that the central cavity is now made 
larger, so that eventually it has the same properties as the dielectric as a whole. 
It seems to me that the expression Kirkwood derives is basically the same as (A. 9) but with 
the dipole moment of the central dipole replaced by M(i?,ro), so that 

^= (. + 2)(2e+l) ^'^-'-°' (S-^) 

and M{R, vq) has the meaning of a dipole moment the central cavity has if it is large enough 
to have a dielectric constant equal to the dielectric constant of the medium as a whole. 
Then Kirkwood argues that actually we may just take the central dipole with its first 
solvation shell, calculate the effective dipole moment of that, which means including the 
local structure. 

In fact this entails taking some limits, so he defines the dipole moment in that limit as 

lim M(R,ro) (B,2) 
so that he can write as a relation for the dielectric constant: 

which then immediately yields eq. (8.10). 



C Macroscopic and microscopic susceptibility. 

In this appendix we discuss some work by Bedeaux et. aL [9, 10] which has some additional 
insights to offer on the second problem mentioned in the final section. It must be remarked, 
however, that these only treat non-polar media. 
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The method is based on the solutions of the full Maxwell equations, which, in k space can 
be written as: 

E{k,uj) = K{k,uj) - F{k,u) • P{k,u;) (C.l) 

In this equation Ey{k,uj) is the solution to the Maxwell equations in vacuum, and F(^,a;) 
is the retarded propagator, given by: 



1 



eo(A;2-(a; + iO)Vc2) 



kk-^ 



(C.2) 



1-1 



We put in an iO to indicate that we give the frequency an infinitesimally small positive 
imaginary part, to insure we only get the retarded solution. 



Derivation of the last equality. 

We need the inverse of the tensor kk — 

To that end we multiply both members with kk — a^, to get 

T • [0 - a^] = 1 

Multiplying on the right with respectively kk and 1 — kk yields 

T . [fcfc - a^] . M - T • kk{k^ - a^) - kk 

and 



T . pfc - a^] • (l - kk) - -a^T • (l - kk) = (l - kk) 
Thus we get for the "longitudinal" and "transverse" components of T: 



T.fcA: = 



kk 



and 



rr. 7T\ I - kk 



Adding these two expressions then gives upon some rearrangement: 
which, with = ]^ — uP- jc? , becomes 



T = 



w2(A;2-a;2/c2) 

This indeed leads to the final equality in (C.2). 
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(C.3) 

(C.4) 

(C.5) 
(C.6) 

(C.7) 
(C.8) 

(C.9) 
(C.IO) 



The next step is to give a constitutive relation between the polarization and the electric field. 
The polarization is given by the sum of the molecular dipole moments (strictly speaking the 
expectation values of the molecular dipole operators): 

p{f,t) = j2fi,mf--m) (c.ii) 

The simplest one can think of is at the molecular level the following 

Mu) = eoX^{u^)'m,cj) (C,12) 

where xi^(^) is the molecular susceptibility tensor of particle i. It obviously depends in 
general on the orientation of the molecule, hence the index i. 

A simpUfication of this is studied by Bedeaux et, ai., where the molecular susceptibility is 
taken to be a scalar quantity, independent of the frequency. In that case the eq. (C.ll) can 
be written as 

P{r, t) = eox"" E Sif- m)E{f, t) = tox"" pif, t)E{f, t) (C.13) 

1=1 

where p(r , t) is the microscopic particle density. 

In the k,uj representation this is not diagonal, but rather a convolution integral, but we can 
write it formally as 

P(fc,a;)-6ox"P0^(^,^) (C.14) 

where we introduced the notation x'^'^O to indicate that we have to deal here with an integral 
operator in fc-space. With this definition (C.l) can be written as 

E{k,u) = Ey{k,uj) - eoF{k,uj) ' x""^ & E{k,Lo) (C.15) 
with the formal solution: 

E{k,uj)= [l+6oF(fc,a;).x'^]"'0^.(^,^) (C.16) 
Averaging this expression yields the macroscopic electric field: 

{E{k,uj)) = ^[l + eoF(fc,a;) • x^^]'') O E^.^) (C.17) 
The macroscopic susceptibility x(^) ^) is defined by the relation 

{P{k^)) = eoxiK^) • (Eik.u;)) = e^x{K^) ' ([l + eoF{Kuj) • x°^]"') O E,{k^) (C.18) 
Averaging (C.14) yields, after substitution of (C.17): 

{P{kM) = (eoX°^ O E{k.^)) - (^oX°" 0 [l + eoF{k.oj) • x"^]'') 0 K{luj) (C.19) 
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Comparison of the previous two equations then gives an expression for the macroscopic 
susceptibility in terms of the microscopic, by eliminating the vacuum field. 



We remark here that we can use the following operator identity, written here for an arbitrary 
operator A: 

[l+A]-'^l-A[l + A]'' (C.20) 

to get 

([l + eoFik.u) . x°^]'') = 1 - F(fc,a;) • ^eoX^^ 0 [l + 6oF(^,a;) • (C.21) 
So that in fact we only need to calculate one average. 

Static fields in an equilibrium dielectric. 

The previous derivation was completely general, and since in this set of notes we restrict 
ourselves mainly to dielectrics in equihbrium, we now specialize to the equilibrium dielectric 
in the presence of a static field. This means that the applied fields become constants in time, 
and all averages are now averages over an equilibrium ensemble of polarizable molecules. 
All time/frequency dependence can thus be neglected, and we get for F(fc): 

m = §, (C.22) 

For a collection of spherical molecules with polarizability a, we find for the microscopic 
susceptibility (see also eq. (3.1): 

eoX°''(r) = ap{f) (C.23) 
A constant vacuum field can be represented in ^-space by: 

E,{k) = {27Tf EJ{k) (C.24) 

To calculate inverses of operators, and subsequent averages, we obviously have to make use 
of series expansions. For the first term in such an expansion we get: 

^ r fJk^ ^ ^ ^ ^ ^ _ 

eoX°P0K(fc)-6o J j^x''''{k-k')E,{k') = ap{k)E, (0.25) 

Upon averaging this becomes 

(P(fc)) = 60 O E,{k)) = i27rfapoEJ{k) (C.26) 

where po is the equilibrium particle density. 

Similarly we can expand (C.18) to lowest order, to get: 

(Pik)) = €ox{k) (l - eoF(^) • x"") 0 (C.27) 
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Performing a similar calculation as the one leading to (C.24), we obtain as an intermediate 
result: 

kk 



{P{k)) = eoX(fc)(27r)3 



1 - apo- 



EJ{k) (C.28) 



eoA;2 

Transforming back to f-space we get for (C.24) and (C.26) respectively: 

(P) = apoE, (C.29) 

and 

<^) = ^»^('-|^)^- 

Combining these last two results, and using x = e - 1, we again find the Cla.usius-Mosotti 
relation: 

= ^ (C.31) 

We also note from eq. (C.28) that the internal field is given by the Lorentz field. So there 
must be a relation between the boundary conditions, and the introduction of correlations 
between the particles, which is something we would like to explore further. 

The second order terms. 

The next term in the expansion of (C.19) gives: 

6^X^"OF(fc)-x^^0^.(^) (C.32) 

The problem with these expressions is that we have to be careful with the interpretation. 
The difficulty is caused by the fact that the F-operator is diagonal in fc-space, whereas x*^^ 
is diagonal in r-space. 

Written in fc-space the full expression becomes: 

^" / (S -*-')■ ■/ ^ - ■ (C.33) 

Introduction of (C.21) and (C,22) then gives: 

The next step is averaging. To that end we need to calculate 

{p{k - k')p{k')) = jdfj df' e-'(fe-fc')-e-^'-' {p{f)pif')) (C.35) 

where 

(p(rXf")) = E {Sir- fi}5{f' - fj)) - po6{r- f') + ^ {6{f- 7^)5{f' - fj)) (C.36) 

Introduction of the pair-correlation function p(^)[ll], then leads to: 

{p{r)p{f')) = po6{f- f') + _ f ) (C.37) 
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Remark. 



Translational invariance is maintained here, but only because the externally applied field 
does not exert a force on the polarizable molecules. Of course there are now the interactions 
between the induced dipoles. This does not lead to breaking of the translational invariance 
though. 

In the papers by Bedeaux et.al a cut-off radius is introduced, which is something we do not 
appear to be needing at all, at least not for the first order term. See the following, however. 

Introduction of this result into eq. (C.33) the gives: 

{p{k-P)p{P)) = (27r)^ +p(2)(fc')] s(k) (C.38) 

Now we do get a divergence, from the first term of this expression. This must somehow 
be related to the fact that an induced dipole produces an infinite field at its own position. 
Infinities like this are customarily avoided by introducing cavities, but somehow it seems to 
me that they should not occur at all, so maybe we can just leave out these terms? (see also 
Fixman [1], who appears to get only the i j terms). 

So let's for the moment subtract this infinite self energy term, and proceed to calculate 
the polarization field. The remainder of the calculation is then relatively trivial. For the 
polarization in (C.19) we get (after back-transformation) up to this order: 

(p)=«..4-f/^.<^'(*-)f .f. (C.39) 

A similar calculation using the remark below (C.19) gives for (C.18) 




(C.40) 



which leads to the Clausius-Mosotti type expression: 

6 + 2 - 360 eoPoJ (27r)3^ ^"h^ ) ^^"^'^ 

It is fairly easy to see that due to the operator identity (C.20), we will in fact get a series 
expansion for the Clausius-Mosotti relation involving higher and higher order correlation 
functions. See also Fixman [1], who gives a shghtly different derivation of this result. 



Neglecting correlations. 

If we neglect correlations, we can replace the n-particle correlation functions p^""^ by (po)"" 
[11]. This then leads to the following series for the Clausius-Mosotti expression: 



e + 2 3eo \ 3eo 
We just note here that this is not equal to the Onsager result for high densities. 
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SOLUTIONS 



tween solute and solvent molecules have hardly any effect in the 
formation of a solution in the gaseous state, because the particles 
are so widely separated. - w 

WHY SUBSTANCES DISSOLVE 

Although it is beyond the scope of this book to mention all factors 
involved in the formation of a solution, the ones of prime importance 
in forming liquid solutions are (1) interaction of solute and solventJ o 
form jolvated part icles. (2) reaction of solute and solvent tn fomo a 
new substance that b ecomes solvated,^ and (3) t he tende ncyjjf syg» 
feinsTo attain a maximum disurdei^*' 



FIG. 10-3 



Solvation. Sojmiwn^k^^ 
solute mojecu whose particles are, 

bpnded together. When water^is used as the solvent^ the process is 
also called aquation or hydration. 

We saw in Chap. 5 that many covalent compounds are composed of 
polar molecules. Water is a famiUar example. When a small crystal of 
an ionic substance such as sodium chloride is placed in water, the 
polar molecules orient themselves about the face of the crystal, as 
shown in Fig, 10-3. The attractive force between the water molecules 



chloride 




water 
molecule 



Polar water molecules and ions 
mutually attract each other. 



and the surface ions is great enough to cause the ions to leave their 
fixed positions in the crystal and to move to positions between the 
water molecules. Water molecules are actually loosely bonded to a 
sodium ion or a chloride ion in solution, as shown in Fig. 10-4. Both 
ions are said to be solvated or aquated. 

Solvation is not limited to ionic compounds. Any polar covalent 
solute may interact with a polar solvent. Glucose sugar is a 

H H H H H 
I I I I I 

0 0 0 0 0 

1 I I I I 

H-C-C-C-C-C-C=0 
I I I I I I 
H H H H H H 

covalent compound with polar groups that dissolves in water 
because of solvation. Each hydroxyl group ( — 0 — H) represents a 



